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I.  Introduction 

Over  the  last  two  decades,  there  has  been  great  interest  in  the  development  of  techniques  to 
generate  and  detect  freely  propagating  pulses  of  coherent  sub-millimeter  wave  radiation  using 
ultrashort  optical  pulses  from  mode-locked  lasers  [1,  2].  These  techniques  use  ultrafast  optical 
excitation  to  produce  either  radiative  current  transients  in  photoconducting  antennas  or  nonlinear 
polarization  via  optical  difference  frequency  mixing  in  a  variety  of  nonlinear  media.  Such 
frequency  down-conversion  schemes  make  it  possible  to  access  the  wide  range  of  frequencies 
within  the  baseband  of  an  ultrashort  optical  pulse,  which  can  easily  extend  out  to  several  hundred 
cm'1  (1  THz  =  33.33  cm'1  =  4.1  meV).  Coupled  with  phase-coherent  detection  using  either 
photoconductive  or  electro-optic  sampling,  these  schemes  of  terahertz  (THz)  generation  have 
provided  powerful  tools  for  time-domain  far-infrared  (FIR)  spectroscopy  [3]  and  THz  imaging  [4]. 

In  studies  supported  by  this  research  grant,  we  have  investigated  both  the  fundamental 
aspects  and  spectroscopic  applications  of  optoelectronic  techniques  of  generation  and  detection  of 
sub-picosecond  electromagnetic  pulses.  The  original  research  objectives  pursued  under  this  grant 
included  the  development  of  significant  new  capabilities  and  an  improved  fundamental 
understanding  of  the  physical  mechanisms  involved  in  the  optoelectronic  generation  and  detection 
of  THz  radiation.  The  characteristics  of  such  optically  excited  THz  emitters  with  respect  to  the 
scaling  of  their  peak  THz  intensities  using  higher  optical  fluences  and  enhancements  in  their  THz 
frequency  response  were  also  to  be  explored.  One  of  the  main  issues  addressed  by  these  studies 
was  the  development  of  new  bright  coherent  sources  and  sensitive  detectors  to  enable  linear  and 
nonlinear  spectroscopy  in  the  far-infrared  region  of  the  electromagnetic  spectrum. 

During  the  period  of  time  covered  by  this  report,  our  accomplishments  include  the 
development  of  ultrafast  laser-based  techniques  and  new  devices  for  generating  and  coherently 
detecting  both  free  space  and  guided-wave  THz  radiation.  The  combination  of  an  optoelectronically 
driven  THz  emitter  and  a  synchronously  gated  optoelectronic  detector,  together  with  appropriate 
FIR  optics,  has  been  utilized  to  measure  the  linear  optical  properties  of  a  variety  of  dielectric  and 
semiconductor  materials  with  exceptional  sensitivity  over  a  spectral  range  exceeding  100  cm'1  [1, 
3].  Because  of  the  gated  nature  of  the  phase-coherent  detection  process,  the  strong  thermal 
background,  which  plagues  traditional  measurements  in  the  far-infrared,  is  observationally  absent 
[3].  As  discussed  below,  significant  progress  has  been  made  towards  enhancing  the  performance 
and  expanding  the  capabilities  of  such  a  system  for  time-domain  far-infrared  spectroscopy. 

II.  Key  Accomplishments 

The  major  accomplishments  of  our  research  efforts  for  the  period  of  November  1,  1991 
through  April  30, 1996  are  summarized  below:- 
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•  Development  of  a  new  large-aperture  Si  p-i-n  diode  for  generating  intense  sub-picosecond 
electromagnetic  radiation.  The  frequency  response  of  such  a  THz  emitter  is  not  determined  by 
the  thickness  of  the  intrinsic  layer  but  by  the  transient  dynamics  of  the  photogenerated  carriers 
within  it.  Thus,  the  device  structure  can  be  optimized  for  high  quantum  efficiency  at  the  pump 
wavelength  without  sacrificing  the  THz  frequency  bandwidth. 

•  Study  of  transient  electron  dynamics  in  GaAs  using  optically  induced  THz  radiation  from 
large-aperture  GaAs  p-i-n  diodes.  The  changes  in  the  amplitude,  duration  and  shape  of  the 
radiated  THz  fields  as  a  function  of  the  wavelength  of  the  optical  excitation  and  the  reverse 
bias  reflect  the  velocity  overshoot  property  of  the  transient  carrier  dynamics  in  GaAs. 

•  Observation  of  Gunn  oscillations  by  triggering  an  appropriately  biased  vertical  transferred 
electron  device  (Gunn  diode)  with  femtosecond  optical  pulses. 

•  Investigation  of  the  dynamics  of  photoinjected  carriers  on  sub-picosecond  time  scales  by 
careful  analysis  of  the  radiated  THz  fields  induced  by  sub-band  gap  optical  excitation  of  bare 
semiconductor  surfaces.  For  such  non-resonant  excitation,  the  major  contribution  to  the 
generation  of  THz  pulses  is  from  the  reactive  photocurrent  caused  by  the  polarization  of  virtual 
carriers  by  the  depletion  field  (field-induced  optical  rectification). 

•  A  systematic  study  of  the  physical  mechanism  involved  in  the  generation  of  THz  radiation  by 
femtosecond  optical  excitation  of  bulk  GaAs.  The  experimental  results  reveal  that  both 
displacement  and  transport  currents  contribute  to  THz  radiation  and  their  relative  contributions 
are  highly  sensitive  to  the  dc  electric  field. 

•  Construction  of  a  Ti:sapphire  regenerative  amplifier  (pumped  by  a  Q-switched  frequency- 
doubled  Nd:YLF  laser)  capable  of  producing  1  mJ,  150  fs  pulses  at  1  kHz.  This  system  is 
being  used  for  high  intensity  broadband  and  narrowband  THz  generation. 

•  Demonstration  of  a  novel  technique  of  generating  tunable  narrowband  intense  THz  radiation  in 
free  space  from  large-aperture  photoconducting  antennas.  This  scheme  relies  on  mixing  two 
linearly  chirped  broadband  optical  pulses  with  a  variable  delay  to  produce  a  tunable  quasi- 
sinusoidal  intensity  modulation. 

•  Generation  and  coherent  detection  of  free  space  narrowband  THz  radiation  using 
photoconducting  dipole  antennas  tunable  from  0.06  to  0.8  THz.  The  spectral  selectivity  of  the 
narrowband  emitter  and  the  narrowband  detector  can  be  tuned  independently. 

•  A  detailed  analysis  of  the  chirped  pulse  beating  process  taking  into  account  limitations 
imposed  by  the  cubic  phase  in  the  chirped  pulses.  The  numerical  analysis  agrees  very  well 
with  cross-correlation  measurements  performed  by  optical  second  harmonic  generation. 
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•  Investigation  of  the  saturation  characteristics  of  the  peak  THz  field  emitted  by  planar 
photoconducting  antennas  with  increasing  optical  fluence.  We  have  shown  that  the  chirped 
pulse  mixing  scheme  can  produce  narrowband  THz  radiation  with  the  highest  possible  spectral 
brightness  from  such  structures. 

•  Generation  of  narrowband  THz  radiation  from  type  II  phase-matched  difference  frequency 
mixing  in  LiNb03  using  femtosecond  optical  excitation.  An  FWHM  bandwidth  of  ~10  GHz 
tunable  over  -400  GHz  was  obtained. 

•  Generation  of  broadband  THz  radiation  from  a  poled  polymer.  The  coherence  length  for 
difference  frequency  mixing  in  this  medium  is  a  few  mm  for  frequencies  >15  THz,  i.e.,  >  100 
times  larger  than  that  of  LiNb03  while  the  FIR  absorption  is  ~10  times  lower.  This  implies 
that  these  materials  are  ideal  candidates  for  intense  FIR  emitters. 

•  The  first  demonstration  of  coherent  detection  of  freely  propagating  THz  radiation  by  electro¬ 
optic  sampling  in  a  poled  polymer.  This  technique  has  the  potential  to  significantly  increase  the 
detection  bandwidth  to  >15  THz. 

•  Demonstration  of  an  ultra-high  bandwidth  system  for  coherent  THz  spectroscopy  based  upon 
the  use  of  non-resonant  second-order  nonlinearities,  i.e.,  optical  rectification  and  electro-optic 
sampling  in  <1 10>  ZnTe  crystals. 

•  Development  and  implementation  of  a  non-perturbative  probe  of  transient  electric  fields  in 
high-  speed  silicon  circuits  and  devices  based  on  optical  second  harmonic  generation.  We  have 
recently  demonstrated  the  first  measurement  an  ultrafast  electrical  pulse  propagating  on  a 
coplanar  transmission  line  with  subpicosecond  time  resolution  by  SHG. 

III.  Optoelectronic  Generation  of  THz  Radiation  from  Semiconductor  Devices 

Substantial  progress  has  been  made  in  the  area  of  generating  broadband  THz  radiation  from 
a  variety  of  semiconductor  devices.  Such  optically  induced  THz  radiation  is  a  powerful  tool  for 
obtaining  valuable  information  on  the  sub-picosecond  dynamics  (such  as  transient  mobility)  of 
photogenerated  carriers  under  applied  fields  in  semiconductors,  and  hot  carrier  effects  associated 
with  them  (such  as  transient  velocity  overshoot).  THz  generation  provides  a  new  optoelectronic 
technique  to  characterize  the  surface  and  interface  electronic  properties  of  semiconductors  through 
the  measurement  of  the  amplitude  and  direction  of  the  surface  static  field,  doping  type  and  dopant 
concentration,  and  even  crystal  symmetry.  Our  work  in  this  area  over  the  past  few  years  is 
summarized  in  the  following  sections. 

A.  THz  radiation  from  large-aperture  Si  p-i-n  diodes: 

One  of  the  difficulties  we  encountered  in  our  previous  work  on  generating  intense  THz 
radiation  from  a  large-aperture  planar  photoconducting  antenna  is  the  limited  attainable  bias  field. 
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The  highest  bias  field  is  determined  not  only  by  the  breakdown  voltage  of  the  semiconductor 
material,  but  also  the  electrode  contacts  and  the  environment  under  which  the  experiment  is 
conducted.  Normally,  the  large-aperture  photoconducting  antennas  we  have  used  before  can 
withstand  only  a  moderate  field  (~  2  kV/cm)  which  is  much  lower  than  the  breakdown  field  of  the 
material.  To  achieve  higher  bias  fields  and  consequently,  more  intense  electromagnetic  radiation, 
we  have  initiated  new  research  on  the  use  of  large  aperture  p-i-n  diodes  for  generating  directional 
and  intense  electromagnetic  pulses  with  terahertz  bandwidths.  By  applying  a  reverse  bias,  we  can 
achieve  a  much  larger  bias  field  in  the  intrinsic  region,  which  is  close  to  the  dielectric  breakdown 
field  of  the  semiconductor. 

The  p-i-n  diode  is  one  of  the  most  common  photodetectors  because  the  depletion  region 
thickness  (the  intrinsic  layer)  can  be  tailored  to  optimize  the  quantum  efficiency  and  frequency 
response.  As  a  photodetector,  high  quantum  efficiency  requires  that  the  intrinsic  layer  thickness  of 
p-i-n  diode  is  larger  than  the  photon  absorption  length.  On  the  other  hand,  the  high  speed  response 
requires  a  thin  intrinsic  layer  for  small  transit  times  of  photogenerated  carriers  through  the  this 
layer.  For  a  conventional  silicon  p-i-n  photodiode  with  3  p.m  thick  intrinsic  layer,  the  3-dB 
response  frequency  is  —10  GHz  for  80%  quantum  efficiency  at  615  nm.  Here,  we  demonstrate  that 
large  aperture  silicon  p-i-n  diodes  can  be  used  as  emitters  for  THz  beams.  This  occurs  when  we 
couple  the  signal  out  by  radiating  into  free  space,  rather  than  coupling  to  the  electrodes,  as  is  done 
in  a  conventional  photodiode.  The  generation  of  optically-induced  THz  radiation  from 
semiconductor  p-i-n  structures  is  similar  to  that  which  we  have  previously  reported  from 
semiconductor  depletion  layers  [5],  It  is  found  that  the  amplitude  and  spectral  bandwidth  of  the 
radiated  pulses  strongly  depends  on  the  reverse  bias  of  the  p-i-n  diode.  Unlike  the  frequency 
response  of  a  conventional  p-i-n  photodetector,  the  frequency  bandwidth  of  the  radiated  pulses 
from  the  large  aperture  Si  p-i-n  diode  is  not  determined  by  the  thickness  of  the  intrinsic  layer  but  by 
the  transient  process  for  the  photogenerated  carriers  to  reach  steady  state  inside  the  intrinsic  layer. 
Hence,  the  radiation  frequency  spectrum  can  extend  to  a  few  THz  and  the  thickness  of  the  intrinsic 
layer  can  be  optimized  for  high  quantum  efficiency. 

The  Si  p-i-n  diodes  used  as  terahertz  radiation  emitters  were  grown  in  a  high  vacuum  Rapid 
Thermal  Chemical  Vapor  Deposition  (RTCVD)  system.  The  carrier  concentration  in  the  intrinsic 
layer  was  less  than  lxlO12  cm"3  and  was  below  the  detection  limit  of  the  spreading  resistance 
measurement.  The  ohmic  contacts  were  fabricated  so  that  an  external  bias  could  be  applied  as 
shown  in  fig.  1.  The  electric  potential  on  the  top  layer  (p+  layer)  is  uniform  due  to  its  small  sheet 
resistance  (1000  Q).  However,  this  sheet  resistance  of  top  layer  is  large  enough  so  that  the 
absorption  of  the  THz  radiation  in  this  layer  is  less  than  5%.  Measurements  of  the  p-i-n  diode  C-V 
characteristics  demonstrated  that  the  capacitance  was  not  dependent  on  the  reverse  bias,  indicating 
complete  and  uniform  depletion  of  the  field  inside  the  intrinsic  layer.  The  p-i-n  junction  broke 
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down  at  about  43  V  reverse  bias  (140  kV/cm  electric  field  in  the  intrinsic  layer).  Since  the 
absorption  length  of  Si  at  615  nm  was  about  2  jim  and  the  depletion  length  in  the  heavily  doped 
semiconductor  was  extremely  small,  almost  all  the  photogenerated  carriers  were  driven  by  the 
depletion  field  inside  the  intrinsic  region  to  form  a  transient  current.  Hence,  we  can  expect  high 
quantum  efficiency  at  this  wavelength. 

Fig.  1  schematically  illustrates  the  set  up  of  the  experiment.  The  Si  p-i-n  diode  is 
illuminated  by  femtosecond  optical  pulses  from  a  balanced  colliding-pulse  modelocked  (CPM)  dye 
laser  which  produces  70  fs  optical  pulses  at  615  nm  with  a  100  MHz  repetition  rate  and  25  mW 
average  power.  The  optical  beam  is  split  into  two  beams  by  a  beam  splitter.  The  unfocused  pump 
beam  illuminated  the  Si  p-i-n  diode  at  an  incident  angle  of  60  degrees.  The  probe  beam  passed 
through  a  variable  time-delay  stage  and  was  focused  on  a  synchronously  gated  dipole  detector  to 
measure  the  temporal  waveform  of  the  radiated  signals.  The  detector  is  placed  1.5  cm  away  from 
the  Si  diode  along  the  direction  of  reflected  optical  beam. 


Figure  1 

The  radiated  electric  fields  detected  from  the  Si  p-i-n  diode  under  different  bias  fields  are 
shown  in  fig.  2  (a).  Some  very  interesting  features  can  be  seen  from  this  figure.  First,  as  the 
reverse  bias  increases,  the  peak  value  of  the  radiation  field  increases  dramatically.  The  peak 
radiation  field  with  40  V  reverse  bias  is  more  than  25  times  larger  than  that  without  external  bias 
and  more  than  50  times  larger  than  that  from  the  surface  depletion  field  of  a  bare  Si  wafer. 
Secondly,  as  the  reverse  bias  increases,  the  pulse  width  of  the  radiation  field  decreases  and  the  time 
delay  corresponding  to  the  peak  of  the  radiation  field  also  decreases.  When  the  p-i-n  diode  is  in 
open  circuit,  the  pulse  width  (FWHM)  is  1.2  ps;  when  the  diode  is  reverse  biased  at  40  V,  the 
radiation  pulse  width  decreases  to  0.76  ps. 
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Time  Delay  (ps) 

Figure  2a 

Fig.  2  (b)  shows  the  normalized  spectrum  of  the  radiation  field  under  the  corresponding 
reverse  biases.  The  radiation  spectrum  clearly  extends  to  higher  frequencies  as  the  reverse  bias 
increases,  which  corresponds  to  a  decrease  in  the  radiation  pulse  width. 


Frequency  (THz) 

Figure  2b 

In  fig.  3,  we  show  the  detailed  bias  dependence  of  the  peak  radiation  field.  The  radiation 
field  increases  linearly  with  the  reverse  bias  at  "low"  fields  (<15  kV/cm  in  the  intrinsic  layer), 
becomes  sublinear  as  the  field  increases,  and  saturates  at  the  extremely  high  fields  (~  130  kV/cm). 
This  saturation  feature  is  almost  the  same  as  that  of  the  steady  state  drift  velocity  of  electrons  in 
silicon  [6].  The  peak  amplitudes  of  the  THz  pulses  have  a  constant  value  in  the  forward  bias 
region,  because  a  large  forward  bias  cannot  be  applied  to  the  intrinsic  region  due  to  the  large 
forward  current.  However,  if  we  extrapolate  the  Ep-V  curve  to  the  forward  bias  region,  we  find 
Ep=0  when  Vo=+0.8  V  which  is  approximately  the  same  as  the  built-in  voltage  (0.85  V). 
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Bias  (V) 

Figure  3 

For  comparison,  we  used  a  Si  p-i-n  diode  with  a  1  pm  thick  intrinsic  layer  as  an  emitter 
under  the  same  experimental  conditions.  We  found  that  the  radiation  signals  were  smaller  than 
those  from  the  3  pm  diode  with  the  same  depletion  field.  The  smaller  radiation  signal  can  be 
explained  by  the  low  quantum  efficiency  because  the  thickness  of  intrinsic  layer  is  smaller  than  the 
absorption  length  of  optical  pulse.  However,  the  shape  of  the  radiation  signal  and  the  field 
dependence  of  radiation  amplitude  for  the  1  pm  p-i-n  diode  are  almost  the  same  as  the  3  pm  diode 
under  the  same  depletion  field.  Consequently,  the  frequency  spectrum  of  the  radiation  field  does 
not  change  with  the  thickness  of  the  intrinsic  layer. 

To  explain  the  experimental  results  above,  it  is  necessary  to  understand  the  mechanism  of 
the  radiation  generation  by  the  optically  induced  carriers.  The  radiation  field  near  the  large  aperture 
Si  p-i-n  emitter  surface  can  be  expressed  as  [7]: 


where  J  is  the  transient  current  density  induced  by  the  optical  excitation  and  the  depletion  field,  d  is 
the  thickness  of  the  intrinsic  layer,  n  is  the  index  of  refraction  at  the  radiation  frequency  and  0  is 


the  optical  incident  angle.  The  transient  current  density  J  is  given  by: 

7  ql  f  t/2>l  , 

J(x,t)=  J— —exp  -ax-—  vd(t-f,E)dt' 

~  o  V  T  7  (2) 

where  q  is  the  electric  charge  of  carriers,  Iop  is  the  peak  intensity  of  the  optical  pulse,  Hcqo  is  the 

optical  photon  energy,  a  is  the  optical  absorption  coefficient,  x  is  the  time  duration  of  optical  pulse 
(assumed  to  be  Gaussian),  v^  is  the  transient  drift  velocity  of  photogenerated  carriers  and  E  is  the 


depletion  field  which  is  constant  for  the  p-i-n  diode.  The  transient  behavior  of  photogenerated  hot 
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carriers  under  high  electric  fields  is  complicated  and  typically  requires  Monte-Carlo  simulations  to 
obtain  quantitative  results  [8].  If  the  hot  carriers  generated  by  the  optical  pulse  have  large  excess 
initial  energy  (maximum  of  0.88  eV  for  electrons  in  our  experiment),  they  will  have  very  large 
phonon  scattering  rates  due  to  the  energy  dependence  of  momentum  relaxation  time.  Thus  the 
strong  dissipation  of  energy  of  the  carriers  due  to  phonon  scattering  slows  down  the  rise  time  of 
drift  velocity  vd  to  about  a  picosecond  [6,  8].  On  the  other  hand,  the  optical  pulse  duration  %  is  so 
small  (~70  fs)  that  the  transient  drift  velocity  vd  is  virtually  zero  for  t'<  %  .  So  the  transient  current 
can  be  written  as  J  =  Nqvd(t,  E),  where  N  is  the  number  of  carriers  generated  per  unit  area  inside 

the  intrinsic  region  per  optical  pulse.  Then,  the  radiation  field  near  the  emitter  surface  is: 

E  Nqvd(t,E)sin6 

nc0£0|n  cos  0  H-  Vl  —  n2  sin2  0^  (3) 

This  formula  gives  a  simple  relation  for  the  radiation  field  as  a  function  of  carrier  density, 
incident  angle  of  the  optical  pulse,  and  the  transient  velocity  of  the  carriers.  From  eqn.  (3),  the 
radiation  field  near  the  emitter  surface  is  directly  related  to  the  transient  velocity  of  the 
photogenerated  electrons,  so  the  amplitude  of  the  radiation  field  is  proportional  to  the  steady  state 
velocity  of  the  electrons.  This  is  why  the  radiation  amplitude  has  the  same  saturation  feature  as  the 
steady  state  drift  velocity  of  electrons  in  Si,  which  is  shown  in  fig.  3,  and  the  frequency  spectrum 
of  the  radiation  field  does  not  change  with  the  thickness  of  the  intrinsic  layer  of  p-i-n  diode.  It 
should  be  mentioned  that  the  radiation  signal  is  detected  1.5  cm  away  from  the  Si  p-i-n  diode. 
Diffraction  causes  the  radiation  waveform  to  become  a  pulse  instead  of  a  step  function  having  a 
duration  comparable  to  the  rise  time  of  the  transient  velocity  of  the  carriers.  However,  the 
saturation  properties  are  unchanged.  Also,  the  shortening  of  the  radiation  pulse  and  the  decrease  in 
the  time  delay  for  the  radiation  to  reach  its  peak  value  imply  that  the  rise  time  of  the  transient 
velocity  of  the  carriers  decreases  when  the  reverse  bias  on  the  diode  increases.  We  also  can 
estimate  the  radiation  field  from  equation  (2).  By  using  N~3xl010  cm"2  per  pulse,  vd~107  cm/s, 
n(Si)=3.42,  Q=n/3  and  0.1  cm2  optical  beam  size,  the  estimated  value  of  the  radiation  field  near  the 
emitter  is  ~2  V/cm.  If  we  assume  the  rise  time  of  the  transient  velocity  is  1  ps,  the  radiation  field 
after  1 .5  cm  propagation  is  a  pulse  with  0.7  V/cm  amplitude  and  about  1  ps  duration. 

It  is  interesting  to  make  some  comparisons  between  externally  biased  semiconductor 
devices  (including  large  aperture  antennas  [9]  )and  semiconductor  surfaces  with  internal  depletion 
layer  fields  [5]  for  THz  generation.  First,  compared  to  the  semiconductor  surfaces,  the  large 
aperture  p-i-n  diode  is  a  more  efficient  radiation  source  because  almost  all  of  the  photogenerated 
carriers  can  be  confined  inside  the  depletion  field  region  by  controlling  the  thickness  of  intrinsic 
layer.  Secondly,  the  depletion  region  of  the  p-i-n  diode  provides  a  controllable  uniform  electric 
field  while  the  other  two  techniques  cannot.  Thus,  we  can  study  the  time  dependence  of  the  drift 
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velocity  under  different  field  strengths  by  applying  different  reverse  biases.  Finally,  compared  to 
the  large-aperture  planar  photoconducting  antenna  which  can  only  sustain  bias  fields  of  few  kV/cm 
and  needs  a  high  voltage  power  supply,  the  large  aperture  p-i-n  diode  is  able  to  produce  higher  bias 
fields  (-140  kV/cm)  with  much  lower  bias  voltages  (-40  V).  As  the  maximum  possible  radiation 
field  is  limited  by  the  bias  field  at  high  optical  fluences  [10],  the  large  aperture  p-i-n  diode  is 
potentially  capable  of  generating  higher  radiation  fields  at  higher  optical  fluence. 

In  conclusion,  we  have  generated  subpicosecond  electromagnetic  pulses  from  large 
aperture  Si  p-i-n  diodes  by  femtosecond  optical  excitation.  The  changes  in  amplitude  and  duration 
of  these  pulses  under  different  reverse  biases  reflect  the  dependence  of  transient  behavior  of 
photogenerated  carriers  on  the  electric  field. 

B.  Study  of  transient  carrier  dynamics  in  semiconductors  using  optically  induced 
THz  radiation 

The  emission  of  THz  radiation  from  optically  excited  semiconductors  takes  place  on  a 
subpicosecond  time  scale,  and  the  mechanism  behind  it  is  very  complicated.  While  models  based 
on  either  classical  drift  currents  or  field-induced  nonlinear  optical  processes  have  been  proposed,  a 
complete  theoretical  understanding  is  still  lacking.  Nevertheless,  the  temporal  waveform  of  the 
THz  radiation  clearly  reflects  the  transient  physical  processes  in  the  semiconductor.  Therefore,  by 
carefully  analyzing  the  waveforms  of  the  radiation  generated  under  different  experimental 
conditions,  one  can  learn  more  about  the  carrier  dynamics  on  such  short  time  scales. 

1.  THz  radiation  induced  by  sub-band  gap  femtosecond  optical  excitation  o f 
semiconductor  surfaces:-  It  has  recently  been  independently  proposed  by  Yamanishi  [11]  and 
Chemla  et  al.  [12]  that  in  a  dc-biased  multiple  quantum  well  structure  (MQW),  a  THz  electrical 
transient  can  be  produced  by  screening  of  the  bias  field  by  virtual  electrons  and  holes  induced  by  a 
femtosecond  laser  pulse  tuned  below  the  absorption  edge.  This  does  not  require  the  generation  of 
the  real  carriers,  i.e.  there  is  no  energy  deposited  in  the  semiconductor.  The  term  “virtual  electron- 
hole  pairs”  refers  to  the  valence  charge  fluctuations  driven  coherently  by  the  intense  laser  field. 
Since  there  are  no  real  collisions,  this  process  is  coherent,  and  consequently,  very  fast.  It  occurs 
on  a  time  scale  -  (  Eg  -  Acoop)_1,  where  Eg  is  the  band  gap  and  hi Oop,  the  optical  photon  energy. 
Later,  Yablonovitch  et  al.  [13]  extended  this  idea  to  bulk  semiconductors  and  give  it  the  name 
“virtual  photoconductivity”. 

Our  experiments  on  the  THz  radiation  bulk  GaAs  sample  excited  with  the  photon  energy 
tuned  from  above  to  far  below  the  bandedge  has  given  us  more  physical  insight  into  this  nonlinear 
process  [14].  The  THz  radiation  is  induced  from  the  depletion  field  region  at  the  surface  of  the 
sample.  Since  the  real  current  disappears  with  the  decrease  in  the  real  absorption,  by  examining  the 
THz  radiation  induced  by  different  photon  energies,  we  can  make  some  assessment  of  the  relative 
contributions  from  the  two  models  to  the  generation  of  the  THz  signal. 
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Two  tunable  femtosecond  lasers  have  been  used  in  these  experiments.  One  is  a  tunable 
dual-jet  hybrid  mode-locked  Styryl  9  dye  laser  synchronously  pumped  at  76  MHz  by  a  frequency- 
doubled  Nd:YLF  laser.  The  other  is  a  self-modelocked  Ti:  sapphire  laser  cw-pumped  by  an  Argon 
ion  laser.  In  both  cases,  the  laser  produces  100  fs  pulses  with  the  wavelength  tunable  from  -800 
nm  to  -900  nm.  The  spectrum  of  the  pulse  has  a  full  width  at  half  maximum  (FWHM)  of  -12 
meV.  Fig.  4  shows  the  schematic  of  the  experimental  set  up.  The  bulk  GaAs  sample  is  placed  in  a 
cryogenic  dewar.  A  photoconducting  dipole  detector  is  located  at  the  specular  angle  of  reflection, 
about  2  cm  away  from  the  sample  which  ensures  detection  in  the  far-field. 


Gating 
pulse 

Figure  4 

The  measurement  is  done  in  the  same  way  as  described  in  section  HI.A.  One  portion  of  the 
unfocused  laser  beam  is  used  to  illuminate  the  bulk  GaAs  sample  to  induce  THz  radiation.  A  probe 
beam,  which  passes  through  a  variable  delay,  and  is  focused  down  to  gate  the  100  |im  dipole 
detector.  The  temporal  waveform  of  the  radiation  is  mapped  out  as  a  function  of  the  relative  time 
delay  between  the  two  beams. 

One  of  the  important  issues  in  the  experiment  is  the  width  of  the  photon  energy  spectrum 
relative  to  the  band  edge.  Due  to  the  Urbach  tail,  at  room  temperature,  the  band  edge  is  very  broad, 
and  becomes  sharper  at  lower  temperature  [15].  Therefore  it  is  necessary  to  perform  the 
measurement  at  low  temperature.  Normally,  we  cool  the  sample  to  about  100  K.  In  order  to  check 
temperature  effects,  we  also  performed  the  same  measurement  at  different  temperatures.  The 
experiments  are  carried  out  at  1 10  K  on  an  undoped  semi-insulating  GaAs  wafer  with  a  thickness 
of  400  pm.  The  optical  beam  strikes  the  sample  at  an  angle  of  about  40  degrees  with  respect  to  the 
normal.  The  average  pump  fluence  is  low  (-3  nJ/cm2 ),  so  that  density  dependent  effects  and  space 
charge  effects  are  small. 

Fig.  5  shows  several  of  the  radiated  waveforms.  5e  denotes  the  energy  difference  between 
the  photon  energy  at  the  center  of  the  laser  spectrum  and  the  bandgap.  These  waveforms  show 
some  very  interesting  features.  When  the  photon  energy  is  well  above  the  band  gap,  the  signal  has 
a  very  large  amplitude  with  a  negative  main  peak,  and  the  amplitude  does  not  change  much  with  the 
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photon  energy.  As  the  photon  energy  is  tuned  towards  the  band  edge,  the  amplitude  starts  to 
decrease.  It  drops  very  quickly  down  to  almost  zero,  as  the  photon  energy  passes  through  the  band 
edge.  At  this  point,  a  very  dramatic  change  occurs.  The  radiated  waveform  reverses  its  sign  with 
its  main  peak  changing  from  negative  to  positive.  This  sign  reversal  is  very  dramatic  depending  on 
the  sharpness  of  the  bandedge.  At  110  K,  this  happens  when  5e  is  tuned  by  only  about  3  meV. 
However,  at  200  K,  this  happens  rather  slowly.  In  both  cases,  the  sign  reversal  occurs  when  the 
absorption  is  reduced  by  more  than  an  order  of  magnitude. 


Time  Delay  (ps) 

Figure  5 

After  the  sign  reversal,  as  5e  is  further  decreased,  the  waveform  is  unchanged.  However, 
the  amplitude  slightly  increases  and  then  decreases  again  at  a  much  slower  rate  compared  to  the 
decrease  near  the  band  edge.  Fig.  6  shows  the  amplitude  of  the  radiated  signal  versus  the  detuning 
energy.  The  signal  is  still  apparent  even  when  the  photon  energy  is  40-50  meV  below  the  bandgap. 
To  gain  more  insight  from  the  experimental  results,  first,  consider  the  radiation  from  a  real 
photocurrent.  In  the  far  field  region,  the  temporal  waveform  should  follow  the  negative  first  time 
derivative  of  the  transient  photocurrent.  It  can  be  roughly  estimated  with  an  assumption  about  the 
time  dependence  of  the  transient  current.  In  our  case,  the  depletion  width  is  roughly  several  |im, 
and  the  static  depletion  field  is  estimated  to  be  about  a  few  kV/cm.  Wysin  et  al.  made  a  detailed 
Monte-Carlo  calculation  of  the  transient  photocurrent  in  GaAs  [16].  Their  results  suggest  that  for 
GaAs  (under  a  bias  field  of  ~  5  kV/cm)  illuminated  by  a  femtosecond  laser  pulse  with  a  center 
photon  energy  of  1.5  eV,  the  photocurrent  will  rise  very  quickly,  in  about  200-300  fs,  reach  an 
overshoot  value,  then  gradually  decrease,  and  relax  down  to  the  normal  steady  state  value  on  a 
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picosecond  time  scale.  The  first  time  derivative  of  the  transient  current  will  have  a  sharp  positive 
peak  followed  by  a  small  slow  negative  lobe.  Therefore,  the  main  peak  of  the  radiated  signal  will 
be  negative.  Using  the  numerical  values  from  their  paper,  we  estimate  that  the  amplitude  of  the 
detected  radiation  would  be  on  the  order  of  a  few  V/cm  for  excitation  above  the  band  gap.  This  is 
close  to  what  we  detected  in  the  experiment.  For  excitation  near  or  just  below  the  band  edge,  the 
signal  should  drop  very  quickly,  since  the  optical  absorption  decreases  exponentially  [15].  This  is 
consistent  with  fast  decay  of  the  signal  as  5e  is  tuned  through  the  band  edge  as  shown  in  fig.  6. 


8e=hG)0p-Eg  (meV) 

Figure  6 

When  the  photon  energy  is  tuned  well  below  the  band  gap,  the  real  photocurrent 
diminishes.  However,  as  described  by  Yamanishi,  Chemla  et  al.,  and  Yablonovitch  et  al.  [11-13], 
there  is  a  reactive  photocurrent  caused  by  the  polarization  of  the  virtual  photocarriers,  which  can  be 
described  by  a  third  order  nonlinear  polarization: 

pnl  =  xO)(0;0,(0o[,-mop)EJ[Kp)E;(fflop)Edt  W 

n\ 

where  %  (0;  0,  -coop,  G)op)  is  the  nonlinear  susceptibility  due  to  the  inverse  Franz-Keldysh  effect. 

For  a  bulk  GaAs  sample  assuming  an  abrupt  band  edge,  Yablonovitch  has  obtained  the  following 
simplified  expression: 


%(3)(0;0,-o)  ,G)  )  = 


1  e  acn  f  h 
24tc2  E„m~  l  A 


1.46x10' 


where  A  is  the  detuning  energy:  A=Eg-hooop  in  meV  [13], 

This  nonlinear  polarization  has  the  spatial  and  temporal  dependence  of  the  optical  excitation 
pulse.  Clearly  it  is  enhanced  by  the  near  resonance  excitation.  The  far  field  waveform  induced  by 
PNL  should  follow  the  negative  second  time  derivative  of  PNL.  The  radiated  waveform  would  have 
a  large  main  peak  with  two  small  lobes  of  the  opposite  polarity.  The  polarity  of  the  main  peak  is  the 
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same  of  that  of  %(3)  which  is  positive.  Therefore,  the  main  peak  of  the  radiated  waveform  generated 
by  the  transient  photocurrent  should  have  an  opposite  polarity  to  that  generated  by  the  nonlinear 
polarization  model.  Furthermore,  the  current  contribution  is  linear  with  the  optical  absorption,  and 
decreases  exponentially  as  the  photon  energy  is  tuned  through  the  bandedge.  The  sign  reversal  and 
the  different  rate  of  change  in  signal  amplitude  as  a  function  of  the  photon  energy  in  our 
experiments  clearly  suggest  a  competition  between  transient  photocurrents  and  nonlinear 
polarization  as  the  dominant  mechanism  generating  THz  radiation. 

In  all  cases,  some  real  absorption  probably  existed,  due  to  the  tail  states  and  impurities 
within  the  depletion  region.  However,  at  110  K  (20-40  meV  below  the  band  edge  of  GaAs),  the 
optical  absorption  should  have  been  negligibly  small.  To  check  how  far  the  optical  pulse 
propagates  in  the  crystal,  we  made  a  very  thin  sample  with  a  thickness  of  only  100  |im.  We  expect 
that  as  the  optical  pulse  passes  through  the  sample,  it  should  generate  THz  radiation  both  in  the 
depletion  region  at  the  front  surface  and  the  depletion  region  at  the  back  surface.  The  pump  beam 
was  aligned  at  Brewster’s  angle,  so  that  the  THz  signal  generated  at  the  front  surface  was  not 
reflected  by  the  back  surface. 

Fig.  7  shows  some  of  the  waveforms  we  obtained  from  the  thin  sample.  For  8e=-7  meV 
and  -15  meV,  the  waveform  is  clean  and  no  reflection  by  the  back  surface  is  observed.  However, 
when  8e  is  50  meY  below  the  bandedge,  we  see  a  second  peak  emerging  2  ps  later  with  an 
opposite  polarity  to  the  first  peak.  The  time  delay  between  the  two  peaks  exactly  matches  the  time 
that  the  THz  signal  generated  at  the  back  surface  needs  to  travel  back  to  front  surface.  This  strongly 
suggests  that  at  50  meV  below  the  band  edge,  the  absorption  is  negligible,  as  the  effects  due  to  the 
real  carrier  drift  cannot  explain  the  amplitude  of  the  observed  THz  signal.  The  signal  is  then  mainly 
due  to  optical  rectification. 
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To  further  prove  the  point,  we  also  measured  the  signal  dependence  on  crystal  orientation. 
We  changed  the  optical  polarization  and  then  monitored  the  amplitude  change  of  the  radiation 
signal.  For  above  band  gap  excitation,  the  THz  radiation  did  not  show  any  dependence  except  that 
due  to  Fresnel  losses  at  the  interface.  This  was  double  checked  at  room  temperature  with  the  crystal 
rotated  azimuthally  about  its  surface  normal.  No  orientation  dependence  was  found  in  all  cases. 
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Figure  8 

For  sub-bandgap  excitation,  however,  as  we  change  the  optical  polarization,  we  observe  a 
large  modulation  of  the  THz  radiation  amplitude.  This  follows  from  the  fact  that  S  °c  d-PNL, 
where  S  is  the  radiated  field  strength  and  d  is  the  unit  vector  parallel  to  the  dipole  axis.  This  can  be 
expressed  as: 

s  =  F(Y„e„er)[(xw  +xOT)+(3x,.,,sin;e,-.(^nx-”)  (6) 

(cos(0i  +0r)  +  3cos(9i  -0r)))cos2  \|/r 

where  the  indices  i,  and  r  indicate  the  incident  and  refracted  beam  respectively,  0i(r)  is  the  incident 
(refracted)  angle,  \|/  is  the  angle  between  the  optical  polarization  in  the  surface  plane  and  the  crystal 
orientation  (100),  and  F  is  the  corresponding  Fresnel  coefficient.  Experimental  results  for  8e=+6 
meY  and  8e=- 1 8  meV  are  shown  in  fig.  8. 

This  experiment  has  demonstrated  radiation  generated  by  a  sub-bandgap  excitation.  It  also 
shows  a  transition  in  the  dominant  radiation  contribution  due  to  the  two  models  proposed. 
However,  there  is  one  unknown  factor,  i.e.  the  strength  of  the  depletion  field.  For  undoped  semi- 
insulating  GaAs,  it  is  very  difficult  to  determine  the  exact  magnitude  of  the  depletion  field.  The 
sample  we  used  in  the  above  experiment  is  believed  to  have  a  very  small  depletion  field.  We 
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roughly  estimate  it  to  be  ~  1  kV/cm,  by  comparing  the  signal  strength  generated  from  the  sample 
with  that  from  a  p-i-n  diode. 

2.  Transient  velocity  overshoot  effects  observed  in  THz  radiation  from  GaAs  p-i- 
n  diodes:-  We  have  demonstrated  that  p-i-n  diode  provide  useful  methods  for  studying  the 
transient  behavior  of  carriers  in  photoconducting  materials.  In  the  following  section,  we  report  the 
study  of  transient  electronic  processes  by  optically  inducing  terahertz  radiation  from  a  large 
aperture  GaAs  p-i-n  diode.  The  GaAs  p-i-n  diodes  used  as  terahertz  radiation  emitters  in  this 
experiment  were  grown  by  MBE.  After  a  1  (im  n+  GaAs  buffer  layer  was  grown  on  a  n+  GaAs 
substrate,  a  3000  A  intrinsic  GaAs  (<10^  cm'3)  layer  was  grown  on  top  of  the  buffer  layer.  This 
was  followed  by  a  200  A  p+  doped  contact  layer  with  8000  Q  sheet  resistance.  The  ohmic  contacts 
were  fabricated  so  that  an  external  bias  could  be  applied.  Due  to  the  low  doping  level  of  the 
intrinsic  layer,  we  expect  complete  depletion  of  the  intrinsic  region  and  uniformity  of  the  depletion 
field  inside  the  3000  A  intrinsic  layer.  The  p-i-n  broke  down  at  5.0  V  reverse  bias  (210  kV/cm 
electric  field  in  the  intrinsic  layer). 

The  set  up  of  our  first  experiment  for  a  large  aperture  GaAs  p-i-n  diode  is  identical  to  that 
for  a  Si  p-i-n  described  above  [7].  Two  femtosecond  laser  sources  are  used  to  produce  terahertz 
radiation  and  gating  the  THz  detector.  For  the  first  set  of  measurements,  the  optical  excitation 
source  is  a  CPM  dye  laser  operating  at  620  nm.  The  unfocused  pump  beam  illuminates  the  GaAs 
p-i-n  diode,  at  an  angle  of  incidence  of  60  degrees.  The  other  beam  passes  through  a  variable  time- 
delay  stage  and  is  focused  on  the  detector  to  measure  the  temporal  radiation  signals.  The  detector  is 
a  radiation  damaged  silicon-on-sapphire  dipole  antenna. 
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The  time  resolved  radiated  electric  fields  detected  from  the  GaAs  p-i-n  diode  at  room 
temperature  under  different  biases  are  shown  in  fig.  9.  Some  important  features  can  be  seen  from 
this  figure.  First,  as  the  reverse  bias  increases  from  1  V  to  -4.5  V  (which  corresponds  to  an  electric 
field  of  13  kV/cm  and  200  kV/cm,  respectively),  the  peak  value  of  the  radiation  field  increases  by 
about  75%.  Secondly,  as  the  reverse  bias  increases,  the  pulse  width  of  the  radiation  field  decreases 
from  0.62  ps  to  0.57  ps,  and  the  delay  time  for  the  radiation  field  to  reach  its  peak  point  also 
decreases.  The  shift  of  peak  of  the  radiation  pulse  is  about  0.2  ps  in  this  process.  The  spectrum  of 
the  radiation  field  in  fig.  9  extends  to  higher  frequencies  as  the  reverse  bias  increases,  which 
corresponds  to  a  decrease  in  the  radiation  pulse  width. 

Fig.  10  shows  the  radiation  signals  from  the  Si  p-i-n  and  GaAs  p-i-n  respectively  under  the 
illumination  of  the  CPM  laser  (620  nm).The  depletion  field  in  the  intrinsic  layer  of  both  p-i-n 
diodes  is  about  130  kV/cm.  The  radiation  waveforms  from  both  diodes  are  basically  unipolar,  but 
the  radiation  signal  from  GaAs  is  faster  and  its  waveform  has  a  small  but  significant  dip  after  the 
main  pulse  as  compared  to  that  of  Si. 


Time  delay  (ps) 

Figure  10 

Velocity  overshoot  of  electrons  in  GaAs  is  expected  to  be  stronger  when  the 
photoexcitation  is  closer  to  the  conduction  band  edge  [16].  Thus,  it  is  desirable  to  repeat  the 
experiment  at  lower  excitation  photon  energy.  In  the  second  experiment,  the  optical  excitation 
source  is  a  Ti:Al203  modelocked  laser.  The  experimental  set  up  is  the  same  as  the  one  using  the 
CPM  laser. 

Fig.  11  shows  the  radiation  signal  from  the  Si  p-i-n  and  GaAs  p-i-n  diodes  under  the 
illumination  of  Ti:Al203  modelocked  laser  (810  nm)  with  the  same  experimental  set  up,  where  the 
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average  power  of  the  laser  beam  illuminating  the  diodes  is  300  mW.  The  depletion  field  in  the 
intrinsic  layer  of  both  p-i-n  diodes  is  about  130  kV/cm. 


Time  delay  (ps) 

Figure  1 1 

The  amplitude  of  the  radiation  from  GaAs  is  about  8  times  larger  than  that  from  the  Si  diode.  The 
pulse  width  of  the  radiation  from  the  GaAs  diode  is  0.50  ps  which  is  also  the  limit  of  the  detector 
response.  It  is  very  interesting  that  the  radiation  signal  is  insensitive  to  the  bias  on  the  GaAs  diode. 
Compared  to  the  radiation  signals  generated  by  620  nm  illumination,  the  radiation  signal  from  the 
GaAs  diode  now  is  noticeably  faster  and  the  pulse  shape  is  significantly  changed,  while  the 
waveform  of  the  radiation  from  Si  diode  is  almost  unchanged. 

From  the  transport  current  model,  the  radiation  field  near  the  emitter  surface  is  directly 
related  to  the  transient  velocity  of  the  photogenerated  electrons  [7].  Therefore,  the  hot  electron 
transient  dynamics  in  GaAs  determines  the  characteristics  of  the  radiation  signal  from  the  diode. 
The  transient  behavior  of  hot  electrons  photogenerated  by  short  optical  pulses  in  GaAs  under 
electric  fields  has  been  studied  both  experimentally  [17]  and  by  Monte-Carlo  simulations  [16].  An 
important  feature  predicted  by  Monte-Carlo  simulations  is  the  well  known  velocity  overshoot,  i.e. 
under  high  electric  field  and  low  photon  energy,  the  peak  velocity  of  the  electrons  will  exceed  the 
steady  state  drift  velocity  before  relaxation  by  phonon  emission.  But  the  experimental  observation 
of  velocity  overshoot  is  difficult.  The  results  of  previous  experiments  using  high  optical  fluences 
[17]  on  overshoot  response  on  picosecond  time  scales  has  been  interpreted  in  terms  of  charge- 
separation  effects  [18].  Also,  a  uniform  electric  field  is  difficult  to  achieve  in  experiments  in  which 
case  geometrical  effects  may  be  important.  However,  under  our  experimental  conditions,  the 
photogenerated  carrier  density  is  less  than  1015/cm3,  and  the  radiation  pulse  is  linear  with  respect  to 
the  optical  power,  which  indicates  no  band  filling,  carrier-carrier  interaction,  or  space  charge 
effects  on  the  radiation  process  [19].  For  THz  radiation  from  the  GaAs  p-i-n  diode,  we  believe  that 
velocity  overshoot  of  electrons  explains  the  experimental  results. 
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For  GaAs,  the  most  important  features  of  the  velocity  overshoot  are  the  dramatic  increase 
of  the  peak  transient  velocity  with  the  decrease  of  optical  photon  energy  and  the  transient  waveform 
of  the  drift  velocity.  We  expect  a  change  in  the  radiation  waveform  when  we  reduce  the  photon 
energy  of  the  optical  excitation  pulse  from  2.0  eV  (CPM  dye  laser)  to  1.54  eV  (Ti:Sapphire  laser). 
Under  a  130  kV/cm  bias  field,  the  amplitude  of  the  radiation  field  from  the  GaAs  diode  illuminated 
at  810  nm  is  8  times  larger  than  that  at  620  nm  if  we  use  the  amplitude  of  the  radiation  field  from 
the  Si  p-i-n  diode  as  a  reference.  By  considering  the  optical  absorption  lengths  at  different  optical 
wavelengths,  the  radiation  field  per  electron  from  GaAs  with  1.54  eV  photoexcitation  (810  nm)  is 
about  5  times  larger  than  that  with  2.0  eV  photoexcitation  (620  nm).  Thus  the  peak  value  of  the 
average  transient  velocity  of  GaAs  with  1 .54  eV  photoexcitation  should  be  about  5  times  larger  that 
with  2.0  eV  photoexcitation,  which  is  in  agreement  with  the  numerical  simulations.  Although  the 
simulations  extend  to  only  50  kV/cm  bias  field,  our  comparison  is  still  believed  to  be  valid, 
because  of  the  saturation  behavior  of  velocity  overshoot  under  high  bias  fields  at  both  optical 
wavelengths. 

For  the  case  of  radiation  from  the  GaAs  diode  by  620  nm  (2.0  eV)  excitation,  the 
significant  feature  is  the  dip  after  the  main  pulse.  Since  the  radiation  waveform  from  Si  has  no  such 
feature  and  the  corresponding  transient  velocity  is  a  well  known  monotonically  increasing  function, 
we  attribute  this  dip  after  the  main  pulse  to  a  weak  velocity  overshoot.  As  the  depletion  field 
increases  and  the  overshoot  effect  becomes  stronger,  this  feature  becomes  more  obvious  as  shown 
in  fig.  9.  Velocity  overshoot  in  GaAs  is  usually  a  very  fast  process  (less  than  0.5  ps  for  2.0  eV 
optical  excitation),  so  that  the  corresponding  radiation  signal  is  faster  than  that  from  Si,  which 
lacks  the  velocity  overshoot  feature.  As  the  photoexcitation  energy  decreases,  the  dip  after  the  main 
pulse  due  to  velocity  overshoot  becomes  larger,  and  finally,  the  radiation  waveform  becomes 
bipolar. 

Velocity  overshoot  can  also  explain  the  bias,  dependence  of  the  radiation  pulse.  Unlike  the 
case  of  the  Si  p-i-n  diode  [7],  the  amplitude  of  radiation  signal  from  the  GaAs  p-i-n  diode  does  not 
have  the  same  features  as  the  steady  state  drift  velocity  of  electrons  in  GaAs.  It  is  well  known  that 
the  steady  state  drift  velocity  of  electrons  in  GaAs  decreases  with  increasing  electric  field  above  3.5 
kV/cm  at  room  temperature.  At  2.0  eV  optical  excitation  and  room  temperature,  the  onset  of 
velocity  overshoot  starts  at  about  20  kV/cm  [16],  and  the  peak  transient  velocity  increases  slowly 
with  the  field  which  explains  the  increase  of  the  radiation  signal  with  the  increase  of  the  depletion 
field.  When  the  photon  energy  is  1.56  eV  which  is  slightly  above  the  bandgap  of  GaAs,  the 
overshoot  velocity  will  saturate  at  about  5x10^  cm/s  when  the  electric  field  is  above  20  kV/cm 
which  is  about  the  smallest  field  in  the  GaAs  diode  in  our  experiment.  This  is  why  the  amplitude  of 
the  radiation  signal  is  insensitive  to  the  bias. 
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By  using  the  Si  signal  as  a  reference,  we  can  obtain  a  quantitative  picture  of  the  transient 
velocity  of  GaAs  under  photon  excitation  from  the  THz  radiation.  In  the  frequency  domain,  the 
detected  radiation  signal  from  large  aperture  p-i-n  diode  can  be  expressed  as  a  convolution, 
E(co)=R(co)vd(to),  where  R(co)  is  the  response  function  including  the  effect  of  detector  response 
and  the  propagation  of  the  radiation  and  vd(ff>)  is  the  ensemble  average  of  transient  velocity.  R(co) 


is  the  same  for  different  radiation  sources  if  the  optical  set  up  is  identical.  Now,  it  is  easy  to  see 
that  the  transient  velocity  of  GaAs  can  be  expressed  as: 
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Time  delay  (ps) 

Figure  12 

Therefore,  once  we  have  the  radiation  signal  from  the  Si  and  GaAs  diodes  under  identical 
optical  alignment,  we  can  derive  the  transient  velocity  of  GaAs  from  the  transient  velocity  of  Si. 
From  the  Monte-Carlo  simulations  [8]  and  the  radiation  pulse  width  of  Si,  we  assume  the  transient 
velocity  of  Si  under  an  electric  field  of  130  kV/cm  and  optical  excitation  at  810  nm  is  a  step 
function  with  a  0.5  ps  rise  time.  In  fig.  12,  the  calculated  transient  velocity  of  GaAs  under  an 
electric  field  of  130  kV/cm  and  optical  excitation  at  810  nm  is  shown.  The  velocity  overshoot 
feature  is  very  clear  and  the  transient  time  for  velocity  overshoot  is  0.45  ps. 

In  conclusion,  we  have  generated  THz  subpicosecond  electromagnetic  pulses  from  large 
aperture  GaAs  p-i-n  diodes  by  femtosecond  optical  excitation.  The  changes  in  amplitude,  duration 
and  the  shape  of  these  pulses  under  different  optical  wavelengths  of  excitation  and  different  reverse 
biases  on  the  p-i-n  diode  reflect  the  velocity  overshoot  property  of  the  photocarrier  dynamics. 

C.  Study  of  the  physical  mechanism  of  optically  induced  THz  radiation 

Recently,  two  conflicting  explanations  have  been  proposed  for  the  origin  THz  radiation 
induced  by  femtosecond  optical  excitation  of  surface  depletion  layers  of  semiconductors.  One 
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states  that  the  photoexcited  carriers  are  accelerated  by  the  depletion  field  that  exists  near  the  surface 
producing  a  time-dependent  transport  current  that  radiates  as  a  Hertzian  dipole  [5].  The  other  model 
proposes  that  the  correlated  photogenerated  electron-hole  pairs  are  polarized  by  the  dc  electric  field 
and  give  rise  to  a  time-dependent  polarization  in  the  depletion  layer,  whose  time  derivative  is  the 
displacement  current  that  can  also  emit  electromagnetic  radiation  [19].  Yet  another  source  of  the 
displacement  current  exists  in  non-centrosymmetric  crystals  such  as  GaAs  where  a  direct  second- 
order  nonlinear  optical  difference-frequency  mixing  can  also  result  in  "optical  rectification"  of  the 
excitation  pulse  even  in  the  absence  of  the  symmetry-breaking  depletion  field  [20]. 

In  an  MQW  structure,  since  the  carrier  transport  is  prohibited  by  quantum  confinement, 
THz  radiation  induced  should  be  caused  only  by  the  displacement  current  [21],  In  bulk  material, 
however,  the  physical  picture  becomes  more  complicated  since  both  the  transport  and  the 
displacement  currents  exist  and  contribute  to  THz  radiation.  Distinguishing  the  effects  of  the  two 
mechanisms  in  a  bulk  sample  has  been  hampered  by  the  lack  of  systematic  experiments.  One  of  the 
difficulties  in  analyzing  the  THz  radiation  from  bulk  material  is  the  uncertainty  of  the  surface 
depletion  field.  The  experimental  results  vary  for  samples  with  different  surface  fields  [14,  22].  In 
previous  experiments,  for  some  GaAs  samples,  the  THz  waveform  showed  a  sign  reversal  with 
the  photon  energy  tuned  around  the  bandedge,  while  for  some  other  samples,  no  change  in  the 
THz  waveform  was  observed  [22]. 

To  gain  a  better  understanding  of  the  physical  mechanisms  involved,  a  systematic  study  of 
the  dependence  of  THz  radiation  on  the  dc  electric  field  and  the  excitation  photon  energy  is 
necessary.  Here,  we  present  the  results  of  our  recent  experiments  on  a  bulk  GaAs  sample.  The 
MBE-grown  sample  we  have  used  is  shown  as  an  inset  in  fig.  13.  It  consists  of  a  1.2  |xm  thick 
intrinsic  GaAs  layer  sandwiched  between  two  undoped  3500  A  thick  Alo.4Gao.6As  layers,  on  top 
of  which  there  is  a  200  A  GaAs  capping  layer  followed  by  a  50  A  Cr  film.  The  external  bias  field 
can  be  applied  between  the  semi-transparent  Cr  contact  on  the  sample  surface  and  the  n+  GaAs 
(100)  substrate.  A  similar  design  was  used  in  ref.  6  for  study  of  the  field  dependence  of  THz 
radiation  in  MQW  structures. 

The  sample  is  placed  in  a  cryostat,  and  is  illuminated  by  an  unfocused  laser  beam  from  a 
self-modelocked  Ti:  Sapphire  laser  with  a  pulse  duration  of  150  fs,  a  spectral  FWHM  of  7  nm,  and 
a  repetition  rate  of  76  MHz.  The  laser  center  wavelength  is  tunable  around  800  nm,  which  enables 
the  study  of  THz  radiation  at  different  photon  energies.  The  photocarrier  density  is  held  at  around 
1012/cm3  to  avoid  space  charge  effects  in  the  GaAs  layer.  The  induced  THz  radiation  is  measured 
by  a  photoconductive  dipole  antenna  detector  [23,  24],  which  is  optically  gated  by  a  split-off 
portion  of  the  laser  beam  which  passes  through  a  variable  time  delay.  The  detector  is  placed  in  the 
direction  of  specular  reflection,  4  cm  from  the  sample,  thus  ensuring  detection  of  the  far-field 
electromagnetic  radiation.  The  ion-implanted  Si-on-sapphire  detector  had  a  time  resolution  of  0.5 
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ps  and  a  corresponding  frequency  response  extending  beyond  2  THz.  To  increase  the  sensitivity  of 
the  detection,  a  hemispherical  sapphire  lens  was  attached  to  the  sapphire  substrate  of  the  dipole 
detector  [24].  The  radiation  waveform  was  detected  by  recording  the  photocurrent  from  the  dipole 
antenna  while  changing  the  relative  time  delay  between  the  pump  and  the  gating  pulses.  Since  the 
detection  scheme  is  coherent;  both  the  amplitude  and  phase  information  of  the  radiation  are 
recorded.  The  detected  waveform  is  a  convolution  in  the  time  domain  of  the  radiated  field  and  the 
detector  response,  therefore,  changes  in  the  shape  and  the  amplitude  of  the  temporal  waveform 
should  reflect  the  physical  properties  of  the  radiative  source. 

In  the  first  set  of  experiments,  we  fixed  the  laser  center  wavelength  at  785  nm  which  was 
approximately  80  meV  above  the  bandedge  of  GaAs  at  a  temperature  of  120  K.  We  then  varied  the 
bias  field,  and  measured  the  induced  THz  radiation.  The  waveforms  showed  no  noticeable  changes 
in  shape.  Fig.  13  shows  the  amplitude  of  the  THz  signal  versus  the  dc  bias  field.  The  amplitude  of 
the  signal  demonstrated  a  linear  field  dependence  in  the  low  field  region.  When  the  bias  voltage 
was  -0.65V,  which  was  close  to  the  flat  band  condition,  the  THz  signal  almost  disappeared.  This 
indicates  the  THz  radiation  observed  here  is  mostly  caused  by  field-induced  effects,  and  the 
intrinsic  bulk  %(2)  effects  are  small.  In  the  high  field  region,  this  dependence  becomes  sublinear, 
and  the  signal  amplitude  saturates  when  the  field  reaches  about  10  kV/cm.  This  type  of  saturation 
has  been  observed  previously  in  a  GaAs  p-i-n  sample  [25]  and  a  MQW  sample  [21]  and  was 
interpreted  in  terms  of  the  transport  and  displacement  current  models,  respectively. 


-2.0  -1.0  0.0  1.0 
Bias  Voltage  (V) 

Figure  13 


In  the  second  set  of  experiments,  we  fixed  the  bias  field  at  different  values,  and  scanned  the 
laser  wavelength  about  the  bandedge  of  GaAs.  In  fig.  14,  the  solid  curves  are  the  THz  waveforms 
measured  under  a  dc  field  of  about  3  kV/cm  and  a  sample  temperature  of  120  K,  where  be  is  the 
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detuning  energy,  i.e.  the  energy  difference  between  the  photon  energy  and  the  band  gap.  The 
waveform  changes  drastically  with  the  detuning  energy.  When  the  laser  energy  is  tuned  below  the 
band  edge,  the  radiated  waveform  has  a  negative  main  peak.  The  strength  of  the  signal  changes 
with  the  detuning  energy,  and  reaches  a  maximum  around  bandedge.  As  8e  is  increased  further, 
the  negative  main  peak  decreases,  and  a  positive  peak  starts  to  develop.  Interestingly,  when  8e  is 
close  to  one  LO  phonon  energy,  which  is  36  meV  for  GaAs,  the  THz  waveform  has  a  very  small 
amplitude  with  some  fine  structure.  Eventually,  the  waveform  changes  polarity  as  the  positive  peak 
develops  and  becomes  dominant. 


Figure  14 

This  polarity  change  of  the  radiation  waveform  with  the  detuning  energy  is  very  similar  to 
the  case  we  observed  previously  in  an  unbiased  GaAs  wafer  [14].  According  to  the  transport 
current  model,  the  far  field  radiation  should  follow  the  first  time  derivative  of  the  photocurrent.  In 
our  case  the  photocurrent  is  expected  from  a  Monte-Carlo  simulation  to  increase  rapidly  reaching 
an  overshoot  value  in  a  few  hundred  femtoseconds,  and  then  relax  down  to  a  steady  state  value  in  a 
few  picoseconds  [16].  The  degree  of  overshoot  depends  on  the  photon  energy  and  the  dc  electric 
field.  The  radiation  waveform  is  expected  to  change  slightly  with  the  photon  energy  due  to  the 
different  amounts  of  velocity  overshoot.  However,  in  general,  the  first  time  derivative  of  the 
photocurrent  would  have  a  positive  main  peak  followed  by  a  small  negative  lobe.  The  velocity 
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overshoot  will  only  affect  the  small  negative  lobes  in  the  radiation  waveform.  In  experiments,  this 
change  is  less  obvious  because  of  the  relatively  slow  detector  response.  No  polarity  change  with 
respect  to  the  detuning  energy  is  expected  from  this  model. 

The  displacement  current  model  requires  the  radiation  waveform  to  follow  the  second  time 
derivative  of  the  nonlinear  polarization  [19].  According  to  ref.  19,  the  time  dependence  of  the 
nonlinear  polarization  should  have  the  time  envelope  of  the  optical  pulse  which  is  independent  of 
the  photon  energy.  The  second  time  derivative  of  the  polarization  would  give  a  symmetric  temporal 
waveform  with  a  negative  main  peak  and  two  small  positive  lobes.  Therefore,  the  two  models 
predict  two  different  time  dependencies  of  the  radiation  field  and  the  polarity  of  the  radiated 
waveform  has  an  opposite  sign.  Neither  the  transport  current  model,  nor  the  displacement  current 
model  alone  can  explain  such  a  waveform  change.  The  sign  reversal  of  the  radiation  signal 
suggests  a  competition  between  the  displacement  and  the  transport  currents  for  inducing  radiation. 

The  theory  of  the  carrier  transport  during  and  at  early  times  after  the  optical  excitation 
should  be  able  to  describe  the  generation  of  polarized  pairs  as  well  as  the  translational  motion  of  the 
carriers  in  the  dc  field.  Neither  the  traditional  semi-classical  approach  to  carrier  transport  nor  the 
steady-state  quantum  mechanical  treatment  proposed  in  ref.  19  are  capable  of  this.  As  it  has  been 
demonstrated  recently  [26],  one  has  to  perform  a  full  time-dependent  quantum  mechanical 
calculation  to  be  able  to  describe  both  effects  at  the  same  time.  According  to  ref.  26,  the  optical 
excitation  creates  carriers  in  localized  wavepacket  states  rather  than  stationary  eigenstates.  These 
wavepackets  move  in  the  dc  field  and  produce  the  transport  current,  but  they  are  also  polarized  by 
the  field  which  leads  to  the  displacement  contribution  to  the  total  current. 

This  theory  is  qualitatively  consistent  with  our  experimental  findings  because  it  predicts  that 
the  relative  contributions  of  transport  and  displacement  current  effects  depend  on  the  detuning,  the 
dc  field  and  the  pulse  duration.  At  high  positive  detunings  the  transport  part  of  the  current  is  found 
to  be  dominant,  while  at  negative  detunings  there  is  only  the  displacement  current  carried  by  the 
virtual  carriers.  The  theory  predicts  a  crossover  from  displacement-dominated  to  transport- 
dominated  behavior  as  the  excitation  frequency  increases,  which  we  believe  is  the  main  reason  for 
the  observed  sequence  of  waveform  changes. 

In  fig.  14,  the  dotted  curves  are  the  theoretical  fit  to  our  experimental  data.  However,  the 
direct  comparison  of  the  theory  with  the  experimental  data  is  complicated  by  a  number  of  factors. 
The  response  function  of  the  detector  is  not  known,  and  for  the  purposes  of  this  study  we  have 
assumed  it  to  have  the  form: 


R(t)  =  cos 


^7tt 

J! 


exp^-(t/0.7)2],t  <  0 
exp|-(t/1.0)2j,t  >  0 


(8) 
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where  t  is  the  time  in  units  of  ps.  This  response  function  is  an  empirical  fit  to  the  waveform  at  high 
positive  detunings.  Our  choice  of  this  function  is  based  on  the  fact  that  the  response  function 
should  be  the  detected  waveform  radiated  from  a  step-function  photocurrent,  and  the  theoretical 
prediction  that  at  high  detuning  the  transport  current,  which  can  be  approximated  as  a  step 
function,  is  the  dominant  radiation  source.  A  different  choice  of  the  response  function  will  only 
affect  the  details  of  the  theoretical  fit.  The  main  features  such  as  the  polarity  reversal  are  basically 
independent  of  the  response  function.  Since  the  theory  only  treats  free  electrons  [26],  we  included 
scattering  phenomenologically  by  introducing  an  exponential  decay  of  the  signal  with  the  time 
constant  of  150  fs. 

The  theory  also  does  not  include  excitonic  effects  which  can  be  expected  to  enhance  the 
density  of  states  near  the  band  gap.  To  account  for  this  effect,  we  have  introduced  a 
phenomenological  detuning-dependent  enhancement  factor  (see  inset  in  fig.  14)  to  match  the 
amplitude  of  the  experimental  and  theoretical  waveforms.  This  enhancement  factor  appears  to  have 
a  dip  at  approximately  one  LO-phonon  energy  above  the  band  gap  which  roughly  accounts  for  the 
effect  of  scattering  which  is  also  not  included  into  the  model  explicitly  [26].  A  comparison  of  the 
solid  and  the  dotted  curves  in  fig.  14  shows  that  the  theory  reproduces  the  observed  sequence  of 
waveforms  reasonably  well  which  confirms  our  initial  suggestion  that  the  change  in  the  waveforms 
is  due  to  the  competition  between  the  transport  and  the  displacement  current  effects.  There  is, 
however,  a  considerable  discrepancy  between  the  theory  and  experiment  around  one  LO-phonon 
energy  above  the  gap  which  indicates  that  the  onset  of  scattering  requires  a  more  accurate 
treatment. 

One  of  the  important  features  of  the  theory  is  its  claim  that  the  displacement  contribution 
decreases  with  the  dc  electric  field.  The  dc  field  not  only  pulls  the  virtual  carriers  apart,  but  also 
restructures  the  electronic  wavefunctions.  It  turns  out  that  the  net  result  is  the  decrease  of  the 
displacement  current  with  the  dc  field  [26].  On  the  contrary,  the  transport  contribution  increases 
linearly  with  the  field.  Considering  the  scattering  processes  involved,  the  transport  contribution 
may  increase  sublinearly  and  eventually  saturate  in  the  high  dc  field  region  [16].  Because  of  this 
unexpected  non-perturbative  result,  at  high  fields  the  contribution  from  the  transport  current 
dominates.  The  time  dependence  of  the  transport  photocurrent  in  this  case,  according  to  the  Monte- 
Carlo  simulation  by  Wysin  et.  al.  is  not  expected  to  have  dramatic  changes  with  the  photon  energy 
[16].  Therefore,  no  significant  waveform  change  in  THz  radiation  with  the  photon  energy  should 
be  observed.  We  have  also  carried  out  the  same  measurement  at  different  dc  bias  fields.  As  we 
increase  the  dc  field,  the  change  in  the  THz  waveform  about  the  band  edge  gradually  diminishes. 
When  the  dc  field  is  approximately  10  kV/cm,  no  noticeable  change  in  shape  of  the  radiated 
waveform  is  observed.  Fig.  15  shows  the  results  for  a  bias  field  of  10  kV/cm  and  sample 
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temperature  of  120  K.  The  amplitude  of  the  signal  again  shows  an  enhancement  near  the  bandedge 
due  to  excitonic  effects.  This  is  again  in  agreement  with  theory. 


Figure  15 

In  conclusion,  at  a  low  dc  bias  field  we  have  observed  a  dramatic  waveform  change  of  the 
THz  radiation  when  the  photon  energy  is  tuned  around  the  bandedge.  This  waveform  change  can 
be  qualitatively  explained  in  terms  of  a  competition  between  the  radiative  contributions  from  both 
displacement  and  transport  currents  using  a  theory  based  on  a  fully  quantum  mechanical  treatment 
of  the  transient  photoconductive  response  of  GaAs.  At  a  high  bias  field,  the  transport  current 
dominates,  and  no  change  of  the  THz  waveform  is  observed. 

D.  Observation  of  Gunn  Oscillation  by  triggering  a  vertical  Gunn  diode  with 
femtosecond  optical  pulses 

Optical  sources  with  subpicosecond  time  duration  are  commonplace  in  high  speed 
optoelectronics  laboratories.  The  current  challenge  is  to  develop  a  reliable  means  of  converting 
such  ultrafast  optical  pulses  into  a  comparable  high-speed  electrical  transient.  Recently,  we  have 
demonstrated  a  new  method  which  will  enable  us  to  generate  intense  bursts  of  electrical  oscillations 
in  the  THz  region  of  the  electromagnetic  spectrum.  This  method  is  based  on  the  idea  of  triggering 
Gunn  oscillations  [27]  in  a  vertical  Gunn  diode  with  a  femtosecond  optical  pulse.  To  do  so,  we 
designed  a  vertical  transferred  electron  device  (TED),  e.g.  a  Gunn  diode,  with  an  intrinsic  buffer 
layer  between  the  active  region  and  the  cathode.  The  speed  of  a  Gunn  diode  itself  is  normally 
determined  by  the  length  of  the  active  region,  and  the  upper  limit  is  intrinsically  set  by  the  time 
carriers  need  to  transfer  from  the  L  valley  to  the  T  valley  which  has  been  measured  to  be  about  2 


ps. 
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Figure  16 

Initially,  the  vertical  device  is  biased  in  such  a  way  that  the  field  in  the  active  region  is 
slightly  below  the  threshold  for  Gunn  oscillations  to  start.  By  choosing  an  appropriate 
femtosecond  optical  excitation  source,  we  can  inject  a  large  number  of  photocarriers  within  the 
buffer  layer.  Therefore,  the  field  of  the  active  region  is  quickly  switched  above  the  threshold  field, 
causing  a  burst  of  electrical  oscillation.  The  duration  of  the  oscillation  is  determined  by  the  time 
required  by  photocarriers  to  sweep  the  buffer  region,  the  carrier  recombination  time,  etc.  The 
coherent  detection  of  this  electrical  burst  can  be  achieved  by  measuring  the  radiated  field  in  free 
space  using  a  photoconductive  dipole  antenna. 

Fig.  16  shows  our  preliminary  results.  The  detected  signal  is  composed  of  two 
contributions,  one  is  the  broadband  signal  caused  by  the  fast  initial  rise  of  the  photocurrent  in  the 
buffer  region,  the  other  is  the  narrow  band  oscillatory  contribution  due  to  the  Gunn  effect  in  the 
active  region.  When  the  bias  field  is  much  higher  than  the  threshold,  no  Gunn  oscillation  is 
observed  as  shown  in  fig.  16a,  while  in  the  case  where  the  field  is  close  to  the  threshold  voltage, 
an  oscillation  is  observed  as  shown  in  fig.  16b.  In  the  experiment  described  above,  the  Gunn  diode 
could  be  replaced  with  a  resonant  tunneling  diode  (RTD),  an  IMP  ATT,  or  any  other  similar  device 
involving  an  intrinsic  high  speed  process. 

IV.  Tunable  Narrowband  THz  Radiation 

In  spite  of  the  abundance  of  optical  transitions  between  energy  states  of  various  systems 
such  as  rotational  and  vibrational  levels  of  molecules,  phonons,  energy  gaps  in  superconductors, 
etc.,  spectroscopy  in  the  sub-millimeter  wave  region  has  traditionally  lagged  behind  its 
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counterparts  in  other  regions  of  the  electromagnetic  spectrum.  This  is  primarily  because  of  the  lack 
of  tunable  bright  sources  and  sensitive  detectors  [28].  This  has  particularly  affected  the  study  of 
nonlinear  and  transient  phenomena  in  this  region.  Optical  difference  frequency  mixing  has  been 
developed  extensively  as  an  efficient  technique  for  generating  and  detecting  tunable  coherent  FIR 
radiation.  The  frequency  content  of  the  nearly  single-cycle  freely  propagating  THz  radiation 
produced  by  conventional  methods  which  use  broadband  optical  excitation  extends  out  to  a  few 
THz.  This  implies  that  other  than  the  inherent  limitations  imposed  by  the  frequency  response  of  the 
combination  of  the  emitting  photomixer  and  detector,  and  the  optical  bandwidth  available,  no  direct 
control  over  the  frequency  content  of  the  THz  radiation  generated  is  possible. 

While  these  broadband  techniques  are  powerful  and  well  suited  for  many  applications, 
certain  important  studies  require  either  tunable  narrowband  THz  radiation  or  detection.  In 
particular,  identification  and  characterization  of  weak  nonlinear  interactions,  such  as  second 
harmonic  generation  in  the  THz  region  will  benefit  greatly  from  spectral  control  of  narrowband 
radiation.  To  this  end,  we  have  demonstrated  a  novel  method  of  both  generating  and  coherently 
detecting  freely  propagating  narrowband  THz  radiation  tunable  out  to  ~1  THz.  We  describe  below, 
the  basic  concepts,  our  experimental  results,  and  applications  of  our  recent  work  on  tunable, 
narrowband  optoelectronic  sources  and  detectors  of  free  space  THz  radiation. 

A.  Chirped  pulse  beating: 

Recently,  we  reported  a  novel  modification  of  the  existing  technique  of  generating 
broadband  free  space  THz  from  large-aperture  photoconducting  antennas  by  incorporating  a  simple 
frequency  tuning  scheme  into  it  [29].  This  scheme  relies  on  the  quasi-sinusoidal  optical  modulation 
obtained  by  beating  two  linearly  chirped  broadband  optical  pulses  with  a  variable  delay.  The  basic 
scheme  for  producing  narrowband  optical  modulation  at  tunable  THz  frequencies  is  illustrated  in 
fig.  17.  A  broadband  transform-limited  optical  pulse  is  temporally  stretched  using  a  pair  of  parallel 
diffraction  gratings  so  that  its  instantaneous  carrier  frequency  is  linearly  swept  in  time  across  the 
output  pulse  duration.  The  chirped  pulse  is  then  split  into  two  pulses  by  a  50/50  beamsplitter  and 
one  pulse  is  delayed  by  a  variable  interval  x  with  respect  to  the  other  using  a  Michelson 
interferometer  arrangement.  The  two  collinearly  propagating  pulses  are  then  recombined  in  the 
optical  mixer  so  that  their  spectra  overlap  in  the  time-domain.  If  these  pulses  are  linearly  chirped, 
then  at  each  point  in  time  their  frequency  content  differs  by  a  constant  beat  frequency  fo,  which  is 
directly  proportional  to  the  time  delay  x  and  the  amount  of  frequency  chirp  in  the  input  pulse.  Any 
square-law  photomixer  such  as  a  photoconducting  antenna  or  a  nonlinear  x(2)  crystal  can  then 
convert  this  quasi-sinusoidal  optical  modulation  into  freely  propagating  pulses  of  narrowband  THz 
radiation  with  a  center  frequency  fo,  which  can  then  be  tuned  simply  by  varying  the  delay  x. 
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Figure  17 

To  illustrate  the  basic  concept  involved,  we  consider  an  optical  pulse  with  a  spectrum 
centered  at  wo,  with  an  amplitude  A(t)  and  phase  %(t),  whose  electric  field  is  given  by 
Ein(t)  =  A(t)exp(-i%(t)-ico0t).  When  such  a  pulse  propagates  through  a  linear  dispersive 

medium,  it  undergoes  a  frequency  dependent  phase  modulation  given  by 

O(m)  =  <F(co0)  +  x0(co  -  m0)  +  +  p(co  -  ©0)3+...  (9) 

Z[L 

Here,  To  is  the  group  delay  at  coo,  p  is  the  carrier  frequency  sweep  rate,  and  P  is  the  cubic  phase 
modulation  that  becomes  important  in  certain  cases  as  discussed  below.  Our  experimental  scheme 
requires  us  to  be  able  to  introduce  a  controlled  amount  of  linear  chirp  while  keeping  the  quadratic 
(i.e.  p)  and  higher  order  chirp  sufficiently  low.  The  phase  modulated  output  pulse  is  then  given  by 

V  oo  oo 

Eout(t)  =  —  J  d(oe_i<0*ei<I,(“)  J  dt'ei(“-“o),A(t')eiz(,')  (10) 

2tt 


In  order  to  be  able  to  perform  analytical  calculations,  we  consider  a  transform-limited  optical  pulse 
with  a  Gaussian  envelope  of  e_1  half  width  o  and  a  Gaussian  spectrum  centered  at  coo  (rad/s),  i.e., 
Ein(t)  =  E0exp(-t2/a2  -  ico0t).  The  output  of  the  Michelson  interferometer  through  which  the 

stretched  Gaussian  pulse  is  passed  is  then  given  by  Itotal  (t)  =  |E0Ut  (t  + 1/2)  +  Eout  (t  -  t/2)|2  where  x 


is  the  delay  between  the  two  arms.  The  output  intensity  can  be  written  as 

1,-W  =  i*(t) + r(t) + 2e{  — V’ ' cos[  + o> J  (id 

where  Gn  =g( l  +  4/p2G4)'/2  is  the  stretched  pulsewidth  and  Gn»G.  The  first  two  terms  in  the 


above  expression  are  the  low  frequency  (dc)  components  of  the  total  intensity  given  by 
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I^t)  =  Eo(cr/on)exp[-2(t ±x/2)2/a^j.  They  correspond  to  the  slowly  varying  Gaussian 

envelope  of  Itotai(t)-  The  last  term  is  the  cross  term  which  contains  the  quasi-sinusoidal  optical 
modulation  at  the  beat  frequency  fo  (Hz)  which  is  given  by 

f.M  =  — ejp  (12) 

7tGnG  27C 

As  can  be  seen  from  the  Fourier  transform  of  WaiCO*  the  bandwidth  of  the  optical  intensity 
modulation  at  fo  is  also  proportional  to  the  linear  chirp  rate  p.  It  is  given  by  the  optical  bandwidth 
Avop,  divided  by  the  pulse  stretching  ratio,  i.e., 


V8  __  V2Gp 


=  Avw  — 


Also,  it  should  be  noted  that  the  cross  term  falls  as  exp(-T2/2G2)  with  increasing  the  delay 
t,  and  this  determines  the  tuning  range  of  the  beat  frequency  fo  to  be^/l/an ,  which  corresponds  to 
the  optical  bandwidth  of  the  input  pulse. 

There  exist  several  techniques  of  introducing  controlled  amounts  of  linear  frequency  chirp 
in  broadband  optical  pulses  using  angular  dispersion  of  linear  optical  elements  such  as  diffraction 
gratings,  prisms,  etc..  A  compact  table-top  setup  consisting  of  a  pair  of  parallel  diffraction  gratings 
can  be  used  to  stretch  -100  fs  optical  pulses  by  factors  ranging  from  10  to  20,000.  The  negative 
group  velocity  dispersion  (GVD)  introduced  by  such  a  grating  pair  is  given  by  [30] 

.  .i  ,  f  r/  x  n2'|-3/2 


1  _  92<F(g))  _  -47t2cb 


-smy 


where  d  is  the  grating  constant,  g  is  the  angle  of  incidence  of  the  input  beam,  and  b  is  the 
perpendicular  separation  between  the  gratings  which  can  be  adjusted  for  the  required  pulse 
stretching  factor  (fig.  17).  This  simple  low-loss  method  of  introducing  large  linear  chirps  has 
inherent  limitations  in  that  higher  order  dispersion  terms  become  significant  for  large  input 
fractional  bandwidths  and  pulse  stretching  factors.  The  cubic  phase  term  in  eqn.  (9)  for  the  grating 
pair  has  been  shown  to  be  [31] 

p_l<d3$(o))N  _  1  l  +  (27tc/©0d)siny-sin2y  ^ 

H  ^G)3  l-((27tc/co0d)-siny)2  J  (15) 

The  analytical  evaluation  of  the  integral  in  eqn.  (10)  considering  terms  up  to  the  third  order 
in  eqn.  (9)  for  O(0))  for  a  Gaussian  pulse  input  gives  Airy  functions  and  is  impossible  for  arbitrary 
pulse  shapes  and  higher  than  cubic  phase  terms.  We  have  numerically  evaluated  eqn.  (10)  using 
discrete  Fourier  transforms  of  the  optical  electric  field.  The  effect  of  higher  order  dispersion  on  the 
output  of  the  chirped  pulse  beating  process  is  to  broaden  the  spectral  width  of  the  intensity 
modulation  at  the  beat  frequency  fo  and  to  reduce  the  effective  tuning  range  to  a  fraction  of  the 
input  optical  bandwidth. 
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1.  Optical  Cross-correlation  measurements:-  We  have  experimentally  verified  the  optical 
intensity  modulation  produced  by  beating  chirped  pulses  with  a  variable  delay  using  standard 
second  harmonic  generation  (SHG)  based  optical  cross-correlation  techniques  [32].  In  our 
experiments,  we  used  -160  fs  (FWHM)  pulses  at  800  nm  (AAfwhm  ~8.3  nm)  at  a  76  MHz 
repetition  rate  from  a  self-modelocked  Ti:sapphire  laser  as  the  input  to  our  pulse  stretcher.  The 
stretcher  consisted  of  two  parallel  diffraction  gratings  (1800  1/mm)  which  were  double-passed  to 
produce  pulses  with  an  FWHM  ranging  from  20  to  200  ps.  After  passing  through  the  Michelson 
interferometer  (fig.  17)  the  stretched  pulse  was  cross-correlated  with  the  output  of  the  Ti: sapphire 
laser. 


We  performed  background-free  intensity  cross-correlation  measurements  using  non- 
collinear  phase-matched  SHG.  We  used  a  0.4  mm  thick  type  I  KDP  crystal  in  our  experiments. 
Fig.  18(a)  shows  the  results  of  a  typical  cross-correlation  measurement  for  a  delay  t  =  3.0  ps 
between  the  chirped  pulses.  The  center  frequency  f0,  the  FWHM  bandwidth  Af,  and  the  relative 
amplitudes  at  f0  of  the  intensity  modulation  produced  by  chirped  pulse  beating  for  various  delays  x 
were  obtained  from  numerical  Fourier  transforms  of  the  recorded  waveforms.  Fig.  18(b)  shows 
the  variation  of  f0  and  Af  with  the  delay  t  and  fig.  18(c)  shows  that  of  the  relative  amplitudes  of 
the  peak  of  the  intensity  modulation  spectra  at  f0. 


Time  (ps)  Delay  (ps) 

Figure  18(a)  Figure  18(b) 

Numerical  simulations  using  discrete  Fourier  transforms  (DFT)  were  carried  out  to  obtain 
intensity  modulation  waveforms  produced  by  beating  chirped  pulses  with  both  quadratic  and  cubic 
phase  modulation  produced  by  the  diffraction  gratings.  The  slope  of  the  experimental  tuning  curve 
(f0  vs.t)  was  used  to  obtain  a  value  of  the  linear  chirp  factor  [t  =2.72x1 023  rad2/s2  from  eqn.  (5). 
Eqn.  (6)  for  the  quadratic  phase  modulation  of  our  grating  pair  then  gave  an  effective  separation  of 
b  =18.3  cm.  The  measured  value  of  b  =18.8  cm  implies  that  the  input  to  the  grating-based 
dispersive  line  had  a  small  amount  of  positive  GVD  as  can  be  expected  from  propagation  through 
various  optics.  We  measured  an  autocorrelation  width  of  -225  fs  for  the  output  of  the  Ti:sapphire 
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oscillator  while  a  cross-correlation  measurement  of  the  stretched  pulse  gave  a  width  of  ~93  ps.  The 
above  value  of  |X  corresponds  to  a  width  ~90  ps  for  a  linearly  chirped  Gaussian  pulse  with  a  8.3 
nm  FWHM  bandwidth  at  800  nm. 

Cubic  and  higher  order  phase  modulation  implies  a  variation  in  the  instantaneous  difference 
frequency  f0  between  the  two  chirped  pulses  across  the  duration  of  the  chirped  pulse  beating 
output,  Itotai(t).  Although  the  duration  of  the  stretched  pulse  is  not  affected  appreciably,  this  leads 
to  a  broadening  of  the  intensity  modulation  spectrum  which  increases  with  the  delay  %  as  can  be 
seen  from  the  fit  to  the  Af  vs.  t  curve  in  fig.  18(b).  However,  the  variation  of  the  beat  frequency  f0 
with  the  delay  %  remains  unaffected. 


Figure  18(c)  Figure  19(a) 

The  figure  also  shows  the  variation  of  Af  with  T  in  the  absence  of  a  cubic  phase  term  in  the 
grating  dispersion  equation  which  is  in  very  good  agreement  with  eqn.  (5).  Another  deleterious 
effect  of  cubic  phase  modulation  borne  out  by  both  the  experimental  data  and  the  simulation  is  the 
decrease  in  the  effective  tuning  range  of  the  beat  frequency  f0.  As  shown  in  fig.  18(c),  the  value  of 
Tluax  at  f0  falls  off  much  more  rapidly  with  the  delay  x  in  the  presence  of  cubic  phase  and  it  no  longer 
follows  the  exp(-T2/2a2)  trend  as  expected  in  the  zero  cubic  phase  case. 

B.  A  tunable  THz  beam  system  using  photoconducting  dipole  antennas:- 

As  discussed  above,  photoconducting  dipole  antennas  can  act  as  broadband  nonlinear 
photomixers  to  convert  the  narrowband  intensity  modulation  produced  by  chirped  pulse  beating 
into  tunable  narrowband  FIR  radiation.  Since  the  radiated  field  is  directly  proportional  to  the 
instantaneous  optical  intensity,  we  expect  that  the  frequency  content  of  the  THz  produced  to  follow 
that  of  the  incident  optical  intensity  I,ota](t).  The  experimental  setup  used  to  generate  and  coherently 
detect  this  tunable  THz  radiation  shown  in  fig.  19(a)  is  similar  to  the  THz  time-domain 
spectroscopy  system  of  Grischkowsky  et  al..  [3].  It  consists  of  identical  100  pm  receiving  and 
transmitting  antennas  fabricated  on  ion-implanted  silicon-on-sapphire  (SOS).  In  the  first  set  of 
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experiments,  the  photoconductive  gap  of  the  transmitting  antenna  is  gated  by  the  output  of  the 
Michelson  interferometer  used  to  mix  chirped  optical  pulses  (A,0  =800  nm,  AX  -8.3  nm),  as 
described  above.  This  dipole  antenna  emits  narrowband  THz  radiation  into  free  space.  This 
radiated  field  is  then  coherently  detected  after  propagation  through  ~40  cm  of  air  by  an  identical 
receiving  dipole  antenna  whose  active  gap  is  gated  by  the  output  of  the  Ti:sapphire  oscillator. 

Fig.  19(b)  shows  a  typical  time-resolved  measurement  of  the  radiated  field  obtained  from  a 
100  Jim  dipole  generator-detector  combination  where  both  the  THz  emitter  and  the  detector  were 
excited  by  -160  fs  pulses  from  the  Ti:sapphire  oscillator.  Fig.  19(c)  shows  a  similar  measurement 
obtained  for  the  case  of  narrowband  THz  produced  by  chirped  pulse  beating  in  the  emitter.  In  both 
the  cases,  the  average  optical  power  in  the  pump  and  the  probe  beams  and  the  bias  on  the  emitter 
were  identical.  The  impulse  response  of  the  overall  THz  beam  system  is  given  by  the  Fourier 
transform  of  the  broadband  THz  field  of  fig.  19(b).  The  high-frequency  part  of  this  response  is 
limited  primarily  by  the  finite  risetime  of  the  photocurrent  in  the  emitter,  the  antenna  response, 
which  drops  when  the  wavelength  of  the  emitted  radiation  becomes  comparable  to  the  length  of  the 
antenna,  and  the  non-ideal  imaging  properties  of  the  THz  optics  [3].  This  limits  the  overall  tuning 
range  achievable  by  the  present  system  to  -0.06-0.8  THz.  Another  factor  limiting  the  signal-to- 
noise  ratio  in  our  experiments  is  the  fact  that  the  peak  signal  falls  by  a  factor  on  the  order  of  the 
pulse  stretching  ratio  (an/a  -800  in  our  case),  as  can  be  seen  from  figs.  19(b)  and  19(c). 

The  center  frequency  f0  and  the  FWHM  bandwidth  Af  of  the  narrowband  THz  were 
obtained  from  Gaussian  fits  to  the  Fourier  transforms  of  the  waveforms  recorded  using  the  above 
setup  for  various  delays  x  between  the  chirped  pulses  incident  on  the  dipole  emitter. 


delay  (ps)  delay  (ps) 

Figure  19(b)  Figure  19(c) 

Numerical  fits  to  the  observed  variation  of  both  f0  and  Af  with  x  were  calculated  for  the  case  of 
chirped  pulse  beating  in  the  presence  of  cubic  phase  modulation  in  the  dispersive  line.  The  values 
of  the  quadratic  and  cubic  phase  terms  (ji  and  (3)  used  for  the  numerical  fits  to  this  data  were  within 


10%  of  those  used  for  the  SHG  cross-correlation  data  in  fig.  19(b).  In  order  to  ensure  that  our 
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emitter  was  not  being  saturated  due  to  screening  of  the  bias  field  by  the  injected  carriers,  we 
verified  the  linearity  of  the  radiated  field  with  the  incident  optical  intensity  [33]. 

C.  Tunable  THz  detectors:- 

In  the  second  set  of  experiments,  the  frequency  tunability  was  incorporated  into  the 
receiving  antenna  instead  of  the  emitter.  This  was  accomplished  by  using  the  narrowband  intensity 
modulation  produced  by  chirp  pulse  beating  to  modulate  the  conductivity  of  the  detecting  dipole. 
The  emitter  was  pumped  by  femtosecond  pulses  from  the  Ti:sapphire  oscillator  while  the  detector 
was  gated  by  the  output  of  the  Michelson  interferometer  described  above.  The  average 
photocurrent  in  the  detector  is  given  by  the  convolution  over  time  of  the  THz  field  and  its  time- 
dependent  conductivity.  Hence,  the  overall  system  response  at  any  THz  frequency  to  can  be  written 
in  terms  of  the  average  photocurrent  recorded  by  the  detector  as 

J.,»  =  (16) 

where  HpC(a>)  is  the  photocurrent  response  function  of  both  the  emitter  and  the  detector  (which 
includes  the  THz  optical  system  response)  [3].  By  reversing  the  roles  of  Ipufnp(co)  and  Iprobe(co)  in 
the  THz  time-domain  spectroscopy  setup  used  in  the  previous  set  of  experiments,  we  now  have  a 
coherent  narrowband  detector  whose  tunability  is  limited  by  its  photoconductive  response  Hpc(a>). 
It  should  be  noted  that  this  scheme  should  work  equally  well  for  an  asymmetric  system  where  the 
emitter  and  the  detector  have  different  impulse  responses  H(co). 


0  10  20  0  500  1000 


Delay  (ps)  Frequency  (GHz) 

Figure  20(a)  Figure  20(b) 

The  results  of  the  time-resolved  measurements  using  the  tunable  narrowband  detector  are 
identical  to  those  obtained  with  a  broadband  detector  (fig.  19(c)).  Here  again,  as  shown  in  fig. 
20(a),  the  center  frequency  and  bandwidth  of  the  detected  portion  of  the  broadband  THz  spectrum 
radiated  from  the  dipole  emitter  show  the  expected  variation  with  the  delay  x  between  the  two 
chirped  pulses  gating  the  detector.  Fig.  20(b)  shows  the  relative  amplitudes  of  the  narrowband 
THz  spectra  detected  by  our  tunable  detector  at  various  values  of  the  center  frequency  (delay  x), 
along  with  the  amplitude  spectrum  of  the  detected  waveform  (dotted  line).  Since  the  average  probe 
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power  remains  fixed,  we  expect  this  narrowband  detector  frequency  response  to  be  a  product  of 
the  variation  of  1^(4)  shown  in  fig.  18(c)  and  the  overall  THz  system  response  obtained  from  the 
Fourier  transform  of  fig.  19(b).  We  have  obtained  reasonable  agreement  between  our  numerical 
calculations  of  this  modified  THz  system  response  and  our  narrowband  THz  data  both  for  the 
tunable  emitter  and  tunable  detector  (fig.  20(b)). 

Thus,  we  have  a  complete  characterization  of  the  frequency  response  of  our  tunable 
narrowband  THz  beam  system  based  on  chirped  pulse  beating  in  either  the  emitting  or  detecting 
dipole  antennas.  This  system  can  now  be  used  for  narrowband  time-domain  spectroscopy  over  a 
wide  frequency  range  in  the  THz  region.  It  should  be  noted  that  the  frequency  resolution  of  this 
technique  is  the  same  as  that  of  conventional  Fourier  transform  spectroscopy,  as  they  are  both 
based  on  a  scanning  delay  line  whose  resolution  is  determined  by  the  reciprocal  of  the  total  time 
scan.  This  technique  may  be  particularly  appropriate  to  nonlinear  spectroscopy  and  the  time- 
domain  study  of  nonlinear  interactions  in  the  FIR,  since  the  emitted  and  the  detected  frequencies 
need  not  be  the  same. 

A  less  fundamental  limitation  is  the  influence  of  cubic  and  higher  order  phase  terms 
introduced  by  the  diffraction  gratings  into  the  broadband  input  pulse.  It  has  been  demonstrated  that 
a  pair  of  prisms  can  be  used  to  compensate  for  the  cubic  phase  distortion  produced  by  the 
diffraetion  gratings  [34],  so  that  improvements  to  the  phase  characteristics  of  the  stretched  optical 
pulse  are  clearly  possible. 

D.  Generation  of  intense  tunable  THz  from  photoconducting  antennas 

We  have  demonstrated  that  large-aperture  photoconducting  antennas,  where  the  linear 
dimension  of  the  illuminated  area  is  much  larger  than  the  center  wavelength  of  the  radiated  field, 
have  made  it  possible  to  produce  optically  steerable  diffraction-limited  beams  of  THz  radiation. 
Further,  due  to  their  size,  large-aperture  antennas  can  sustain  very  high  optical  pulse  energies  and 
high  bias  fields,  and  are  hence,  ideally  suited  to  scale  up  the  peak  power  of  the  free  space  THz 
pulses  without  loss  of  THz  bandwidth.  The  availability  of  such  high  energy  far-infrared  pulses 
synchronized  with  intense  subpicosecond  optical  pulses  from  table-top  amplifiers  has  enabled 
numerous  time-resolved  optical  pump  and  THz  probe  measurements  such  as  the  study  of  relaxation 
dynamics  of  photoexcited  carriers  in  GaAs  and  InP.  However,  it  is  now  well  known  that  the 
radiated  fields  from  these  broadband  emitters  saturate  with  increasing  optical  fluence  due  to  the 
screening  of  the  bias  field  by  the  injected  carriers  [33,  35]. 

For  all  our  high  optical  fluence  experiments,  we  have  constructed  a  Ti:sapphire  regenerative 
amplifier  which  is  seeded  by  the  self-modelocked  Ti:sapphire  oscillator  described  above  and 
pumped  by  a  frequency-doubled  Q-s witched  Nd:YLF  laser.  This  system,  based  on  the  concept  of 
chirped  pulse  amplification  to  avoid  deleterious  nonlinear  effects  in  the  gain  medium,  is  capable  of 
producing  1  mJ,  150  fs  pulses  (corresponding  to  peak  power  >  5  GW)  at  1  kHz.  We  have 
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explored  the  physics  of  saturation  mechanisms  in  photoconducting  antennas  made  up  of  various 
semiconductors  by  studying  their  behavior  under  such  intense  broadband  and  tunable  narrowband 
optical  excitation.  We  have  recently  shown  that  the  use  of  stretched  optical  pulses  in  conjunction 
with  such  large  aperture  antennas  can  enable  the  achievement  of  higher  optical-to-electrical 
conversion  efficiencies  within  a  given  FIR  bandwidth.  This  arises  from  the  possibility  of  utilizing 
much  higher  optical  fluences  without  deleterious  saturation  effects  due  to  the  bias  field  limited 
boundary  conditions  at  the  surface  of  the  THz  emitter. 

Extensive  numerical  simulations  of  the  saturation  in  planar  photoconducting  antennas  using 
a  simple  model  for  the  radiated  field  from  a  time- varying  photocurrent  [36]  have  shown  that  the 
spectral  brightness  of  the  narrowband  emitters  can  exceed  that  of  broadband  emitters  by  factors  as 
high  as  200.  The  degree  of  enhancement  over  passive  filtering  of  broadband  signal  is  predicted  to 
increase  with  decreasing  carrier  lifetime  (photocurrent  relaxation  time).  The  results  of  the 
experiments  performed  with  both  large  aperture  and  dipole  antennas  have  demonstrated  the 
possibility  of  working  in  this  saturation  regime.  The  observation  that  the  peak  THz  signals 
obtained  for  a  given  optical  pulse  energy  do  not  scale  with  the  pulse  stretching  ratio  as  expected 
from  a  simple  first  order  calculation,  is  a  clear  indication  of  the  advantage  of  using  narrowband 
emitters.  Our  technique  provides  an  approach  uniquely  suited  to  the  generation  of  pulses  of  intense 
narrowband  THz  radiation  tunable  over  a  wide  range  of  frequencies  optically  synchronized  with 
very  intense  broadband  optical  sources.  This  capability  can  be  exploited  to  probe  FIR  excitations  of 
surfaces  and  interfaces  using  nonlinear  spectroscopic  techniques  such  as  sum-frequency 
generation.  In  the  area  of  far-infrared  nonlinear  optics,  various  methods  of  harmonic  generation 
can  be  readily  explored  with  our  technique. 

E.  Tunable  far-infrared  radiation  from  difference-frequency  generation  in  lithium 
niobate:- 

Optical  difference  frequency  generation  (DFG)  using  broadband  optical  excitation  of 
nonlinear  dielectric  crystals  such  as  LiNb03  has  been  demonstrated  as  a  powerful  technique  to 
generate  coherent  free  space  FIR  radiation  [28,  37].  The  frequency  bandwidth  of  the  THz  radiation 
generated  by  optical  rectification  in  such  nonlinear  crystals  is  limited,  in  principle,  only  by  the 
bandwidth  of  the  femtosecond  optical  excitation.  However,  as  discussed  in  the  previous  sections, 
broadband  phase  matching  is  difficult  to  achieve  due  to  the  high  FIR  dispersion  in  LiNb03.  We 
may  turn  this  property  to  an  advantage  by  using  this  frequency-selective  process  of  phase  matching 
[38]  to  control  the  spectral  content  of  the  THz  radiation  produced. 

We  have  used  type  II  phase  matched  difference  frequency  mixing  in  LiNb03  to  generate 
narrowband  THz  radiation  with  -150  fs  pulses  at  800  nm  from  a  self-modelocked  Ti:sapphire 
laser.  The  FIR  refractive  indexes  and  absorption  coefficients  were  measured  by  THz  time-domain 
spectroscopy.  These  measurements  were  applied  to  calculate  frequency-dependent  phase  matching 
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curves  from  0.1  to  1.0  THz.  Using  these  data,  we  have  obtained  FIR  radiation  tunable  from  0.3  to 
0.7  THz  with  a  bandwidth  of  ~12  GHz.  The  measured  THz  bandwidths  agreed  well  with  the 
dispersion-limited  bandwidths  calculated  from  the  experimental  dispersion  data  for  LiNb03.  Fig. 
21  shows  the  observed  THz  field  from  a  5  mm  thick  LiNb03  cut  at  15°  off  the  c-axis  detected 
using  a  100  jim  SOS  dipole  detector.  The  detected  field  spectrum  (fig.  21b)  shows  both  the 
forward  and  the  backward  phase-matched  DFG  signals  generated.  The  same  can  be  seen  in  the 
temporal  waveform  (fig.  21a). 
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Figure  21a  Figure  21b 

We  believe  that  the  application  of  broadband  chirped  pulses  will  enable  us  to  scale  this 
narrowband  THz  field  to  greater  intensities  than  those  possible  with  femtosecond  pulses.  This 
expectation  is  based  on  the  possibility  of  going  up  to  higher  optical  fluences  on  the  nonlinear 
crystal  without  causing  optical  damage  when  stretched  pulses  are  employed.  By  using  nonlinear 
media  such  as  LiNb03  and  organic  crystals  such  as  DAST,  we  also  hope  to  extend  the  tuning  range 
of  our  narrowband  THz  system  to  the  optical  bandwidths  available. 

V.  Ultra-broadband  THz  Radiation 

In  this  section,  we  review  recent  progress  in  our  laboratory  on  the  generation  and  detection 
of  freely  propagating  THz  radiation  by  means  of  second-order  nonlinear  effects.  As  discussed  in 
the  preceding  sections  ,  remarkable  progress  has  been  made  in  the  development  of  spectroscopic 
capabilities  for  coherent  THz  measurements.  Much  of  the  work  reported  to  date  relies  on  the  use  of 
photoconductive  elements  both  as  emitters  and  detectors  [3,  39].  While  these  devices  are  highly 
optimized  in  many  respects,  there  are  fundamental  limitations  in  their  frequency  response 
associated  with  the  natural  time  constants  for  carrier  dynamics.  The  ready  availability  of  extremely 
short  laser  pulses  suggests  the  potential  for  extending  the  bandwidth  of  coherent  spectroscopy  to 
significantly  higher  frequencies.  For  example,  a  10-fs  transform-limited  laser  pulse  at  X0  =  800  nm 
has  a  bandwidth  of  44  THz  or  roughly  1500  cm1.  The  use  of  materials  with  an  instantaneous 
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optical  response,  such  as  the  non-resonant  electronic  second-order  nonlinear  susceptibility  %(2), 
may  permit  this  enormous  bandwidth  to  be  captured  both  for  the  generation  and  detection  of 
ultrafast  electromagnetic  transients. 

A.  Generation  of  THz  radiation  from  poled  polymers: - 

In  recent  years,  terahertz  radiation  has  been  generated  from  a  variety  of  nonlinear  optical 
materials  [40-42].  Many  of  these  materials  exhibit  a  high  degree  of  dispersion,  corresponding  to 
small  coherence  lengths  for  difference  frequency  mixing,  and  high  absorption.  This  limits  the  total 
conversion  efficiency.  Polymeric  media  offer  many  attractive  properties  for  far  infrared  generation. 
In  addition  to  being  relatively  low  loss,  low  dispersion  materials,  they  are  readily  processable,  can 
be  poled  to  induce  non-centrosymmetry,  and  may  possess  nonlinear  coefficients  greater  than  that 
of  lithium  niobate.  As  discussed  below,  we  have  generated  THz  radiation  from  an  organic  poled 
polymer.  We  believe  that  these  materials  have  potential  both  for  very  intense  FIR  emission  and  for 
extending  the  spectral  range  well  beyond  that  attainable  using  photoconductors. 

We  used  a  conventional  pump-probe  experimental  setup  for  generating  and  detecting  the 
THz  radiation  [41]  with  a  poled  polymer  used  as  the  emitter.  The  polymer  consisted  of  a 
copolymer  of  4-N-ethyl-N-(2-methacryloxyethyl)amine-4'-nitro-azobenzene  (MAI)  and  methyl 
methacrylate  (MMA)  obtained  from  AlliedSignal  [43].  In  the  experiments,  a  free-standing  16  mm 
thick  film,  with  dipolar  alignment  perpendicular  to  the  surface  was  used  [44].  The  polar  axis  (z')  of 
the  uniaxial  nonlinear  medium  was  determined  by  the  direction  of  the  poling  field.  The  optical 
excitation  pulses  were  generated  by  a  mode-locked  Ti:sapphire  laser  operating  at  800  nm.  An  s-  or 
p-polarized  pump  beam,  with  average  power  up  to  600  mW  and  spot  diameter  of  1  mm  at  the 
sample,  was  used  to  pump  the  polymer  film.  The  generated  THz  radiation  was  focused  onto  a  100 
pm  radiation  damaged  silicon-on-sapphire  dipole  detector,  aligned  to  detect  p-polarized  radiation. 
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We  first  measured  the  refractive  index  and  absorption  of  the  MA1:MMA  polymer  film.  This 

was  accomplished  by  THz  time-domain  spectroscopy  [3]  using  a  bare  <100>  GaAs  wafer  as  the 

emitter.  Fig.  22  shows  the  resulting  refractive  index  and  power  absorption  from  0.2  to  1 .0  THz. 

As  can  be  seen,  the  absorption  is  very  weak  over  the  entire  spectral  range  and  the  refractive  index 

is  essentially  flat.  We  believe  that  the  observed  oscillations  are  artifacts  due  to  diffraction  effects 

caused  by  the  aperture  and  the  finite  Fourier  transforms  taken  in  the  data  analysis.  The  coherence 

length,  £c,  for  difference  frequency  mixing  is  given  by 
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where  c  is  the  speed  of  light,  nt  and  nopt  are  the  THz  and  optical  frequency  refractive  indexes, 
respectively,  0lnfe  is  the  far-infrared  frequency,  Xopt  is  the  optical  wavelength,  and  dnopt/dA,  is  a 
measure  of  the  optical  dispersion.  The  coherence  length  is  >1  mm  for  all  FIR  frequencies  below  10 
THz.  For  comparison,  the  corresponding  coherence  lengths  in  LiNb03  are  nearly  200  times 
smaller.  If  we  assume  uniform  generation  throughout  a  1  mm  thick  polymer  with  an  absorption 
coefficient  of  a  ~  4  cm'1,  only  10%  of  the  far  infrared  radiation  would  be  reabsorbed  at  frequencies 
between  0.2  and  1.0  THz.  Thus,  we  should  be  able  to  effectively  exploit  the  large  coherence 
lengths  of  these  polymeric  materials. 


Delay  [ps] 

Figure  23 

The  low  frequency  electro-optic  coefficient  of  our  16  |im  thick  sample  was  r33  ~  1 1  pm/V 
[43,  44].  The  peak  signal  from  the  polymer  sample,  shown  in  fig.  23,  was  roughly  a  factor  of  4 
smaller  than  the  maximum  signal  obtained  from  a  LiNb03  crystal  with  a  thickness  of  1  mm.  The 
relatively  small  difference  in  the  peak  THz  signals  observed  between  the  two  materials  is  due  to  the 
short  coherence  lengths  and  large  THz  refractive  index  of  LiNb03,  corresponding  to  greater 
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Fresnel  reflection  coefficients.  Note  that  the  intrinsic  nonlinear  response  (r33  =  31  pm/V)  of 
LiNb03  is  comparable  to  that  of  our  polymer. 

In  order  to  verify  that  the  THz  signal  is  indeed  due  to  optical  rectification,  we  measured  the 
p-polarized  THz  radiation  as  a  function  of  incidence  angle  and  pump  polarization.  The  experimental 
data  obtained  from  the  16  |im  thick  sample  is  shown  in  fig.  23  for  p-  and  s-polarized  pump  beams 
for  angles  of  incidence  0  ranging  from  0°  to  60°.  The  radiated  electric  field  is  proportional  go2Pnl, 
where  ©  is  the  far  infrared  frequency  and  PNL  is  the  effective  nonlinear  polarization.  For  an  s- 
polarized  pump  beam  and  a  p-polarized  THz  beam,  the  radiated  field  is  proportional  to  PNL(pss)  ~ 
A[d31sin0m]ExEx\ 


Rotation  Angle  [degrees] 

Figure  24 

Similarly,  for  a  p-polarized  pump  beam  and  a  p-polarized  THz  beam,  the  radiated  field  is 
proportional  to  PNL(pPP)  ~A'[(3d33+3d3,)sin0m  +  (3d31-d33)sin(30m)]  EzEz*/4,  where  A  and  A'  are 
constants  and  0m  (=  sin'[sin0  /  nopt])  is  the  propagation  angle  in  the  polymer  medium.  The 
experimental  data  shown  in  fig.  24  is  corrected  for  the  angularly  dependent  transmission 
coefficients  and  effective  pathlength.  We  simultaneously  fit  to  the  two  nonlinear  polarization 
expressions  assuming  A=A'.  The  resulting  best-fit  yields  a  value  of  d31/d33  =  0.34.  A  simple 
thermodynamic  model  for  poled  polymers  [45]  yields  a  calculated  value  of  d31/d33=  0.30  [44], 

In  future  experiments,  we  hope  to  demonstrate  scaling  of  the  THz  radiation  to  increasingly 
high  intensities.  This  may  be  achieved  using  polymers  with  of  millimeter  thickness,  as  argued 
above.  Use  of  amplified  laser  pulses  should  further  enhance  the  THz  radiation  intensity.  Unlike 
photoconductively  gated  devices,  the  radiated  electric  field  from  polymers  is  expected  to  vary 
linearly  with  the  pump  intensity  up  to  the  damage  threshold  of  the  medium.  We  also  intend  to 
verify  the  capabilities  of  these  materials  to  support  optical  rectification  of  laser  pulses  as  short  as  10 
fs,  with  corresponding  bandwidths  of  I0’s  of  THz. 
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B.  Coherent  detection  of  free  space  THz  radiation  by  electro-optic  sampling: - 

To  date,  the  only  broadly  applicable  technique  to  coherently  detect  THz  radiation  requires 
the  use  of  synchronously  gated  photoconducting  dipole  detectors  [3,  39].  While  these  detectors 
exhibit  excellent  sensitivity,  it  is  difficult  to  calibrate  the  measured  electric  field.  Additionally,  there 
is  a  strong  speed  versus  sensitivity  trade-off  which  is  determined  by  the  photoconductive  response 
and  antenna  dimensions  [3].  While  electro-optic  sampling  [46]  has  been  used  extensively  to  study 
ultrafast  electrical  pulses,  these  applications  have  been  primarily  limited  to  transmission  line  and 
circuit  characterization.  The  significant  advantages  of  this  technique  are  the  capability  of  detection 
that  is  limited  only  by  the  optical  bandwidth  and  the  ability  to  calibrate  the  detected  electric  field. 
The  former  is  especially  true  with  the  use  of  poled  polymers,  where  the  material  response  time  is 
negligible  [43,  47]. 

Electro-optic  sampling  has  recently  been  demonstrated  as  a  technique  for  the  coherent 
detection  of  freely  propagating  sub-millimeter  wave  radiation  [48,  49].  In  our  experiments 
described  below,  the  emitter  and  detector  are  physically  separated  by  a  large  distance.  The  electro¬ 
optic  sampling  is  performed  in  a  poled  polymer  that  is  similar  to  that  discussed  in  the  previous 
section.  We  fabricated  the  electro-optic  sampling  element  (EOSE),  shown  schematically  in  fig.  25, 
on  an  R-plane  sapphire  substrate. 

Si  ball  — ► 

HR  coating  — » 


electrodes 


.sapphire 

-polymer 


Figure  25 

A  high  reflectivity  dielectric  coating  centered  at  800  nm  was  first  evaporated  onto  one  side 
of  the  substrate.  Two  coplanar  5  mm  x  5  mm  x  2000  A  thick  aluminum  pads  separated  by  50  pm, 
used  for  poling,  were  then  photolithographically  defined  onto  the  coating.  A  10  pm  thick  film  of 
the  nonlinear  copolymer  was  then  cast  from  solution,  baked  at  130°  C  for  3  hours  in  air,  and 
poled.  A  hyperhemispherical  silicon  lens  was  attached  to  the  sapphire  substrate  to  focus  the  THz 
radiation  into  the  electrode  gap. 

The  experimental  setup  used  for  generating  and  detecting  THz  radiation  is  shown  in  fig. 
26.  A  mode-locked  Ti:sapphire  laser  operating  at  800  nm  with  energy  of  10  nJ  and  160  fs  pulse 
duration  was  used  to  generate  and  detect  the  transient  sub-millimeter  wave  pulses.  A  p-polarized 
pump  beam,  chopped  at  95  kHz,  was  used  to  drive  a  large  aperture  photoconducting  antenna.  The 
emitter  was  biased  with  voltages  ranging  from  +300  V  to  -300  V  and  oriented  so  that  the  resulting 
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THz  radiation  was  p-polarized.  A  pair  of  off-axis  paraboloidal  mirrors  was  used  to  collimate  and 
focus  the  THz  radiation  into  the  EOSE.  The  total  separation  between  the  emitter  and  detector  was 
60  cm. 


Figure  26 

The  detection  system  employed  a  crossed  polarizer  arrangement  with  differential  detection 
[46].  The  5  mW  probe  beam  was  optically  biased  at  its  quarter  wave  point  by  a  Soleil-Babinet 
compensator  and  focused  into  the  electrode  gap  of  the  EOSE.  The  c-axis  of  the  poled  polymer  was 
parallel  to  the  incoming  p-polarized  THz  radiation.  The  reflected  probe  beam,  split  by  the 
Wollaston  prism,  generated  a  quiescent  current  of  approximately  1  mA  in  each  detector.  The 
differential  signal  was  amplified,  digitized  and  averaged  using  a  computerized  data  acquisition 
system. 


The  temporal  response  of  the  system  with  an  antenna  bias  of  300  V  is  shown  in  fig.  27. 
We  calibrated  the  THz  response  of  the  EOSE  by  applying  a  DC  bias  across  the  50  pm  electrode 
gap.  Such  a  calibration  requires  that  we  neglect  the  difference  in  field  distributions  between  the 
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diffracted  THz  field  and  the  fringing  applied  voltage.  We  also  neglect  the  dispersion  of  the  electro¬ 
optic  coefficients.  The  latter  assumption  reasonable  since  the  optical  nonlinearity  is  electronic  in 
nature,  as  shown  by  comparing  electro-optic  and  second  harmonic  generation  measurements  [47]. 
The  peak  detected  field  was  found  to  be  ~  18  V/cm  with  a  noise  floor  of  2.9  (V/cm)  /VHz . 

The  observed  DC  electro-optic  coefficient  r33  was  determined  to  be  11  pm/V, 
corresponding  to  a  value  of  EK,  the  field  needed  to  cause  a  %  phase  retardation,  of  12  MV/cm.  In 
the  shot  noise  limit,  we  found  that  the  minimum  detectable  field  for  our  detection  system  is 

E„.=E  .pOX'SE).  (18) 


21,  -,/Hz 


Thus,  for  a  quiescent  current  of  1  mA  and  1  Hz  bandwidth,  the  minimum  detectable  field  is  100 
(mV/cm)/  a/Hz  .  Our  present  detection  noise  floor  is  nearly  29  dB  above  the  shot  noise  limit  and 
appears  to  be  dominated  by  excess  laser  noise.  We  are  currently  working  to  reduce  the  system 
noise. 


In  the  far  field,  the  radiated  THz  electric  field  is  expected  to  have  a  bipolar  shape,  such  that 
total  time  integral  is  zero  [39,  49].  The  observed  shape  of  the  detected  waveform  in  fig.  27, 
however,  is  unipolar.  This  effect  is  the  result  of  the  perturbation  of  the  applied  THz  field  induced 
by  the  metal  electrodes  in  the  EOSE.  A  similar  effect  was  observed  in  the  first  description  of 
photoconducting  Hertzian  dipole  detectors  [50].  We  modeled  the  influence  of  the  electrodes  by 
considering  classical  electromagnetic  diffraction  through  a  2-dimensional  conducting  screen  with  a 
narrow  slit  [50].  Such  a  model  is  appropriate  since  the  electrodes  are  significantly  larger  than  the 
wavelengths  of  interest.  To  first  approximation,  the  diffracted  field  (and  corresponding  electrode 
voltage)  is  proportional  to  the  time  integral  (fffECco))  of  the  incident  field.  This  is  equivalent  to  the 
charging  of  the  gap  capacitance  by  the  incident  THz  radiation.  An  approximate  expression  for  the 
voltage  across  the  slit  can  be  written  as  [50] 


(19) 


where  Z0  is  the  characteristic  impedance  of  the  medium  and  Cs  ,  which  has  a  weak  frequency 
dependence,  is  the  static  capacitance  of  the  electrode  structure  per  unit  length. 

We  demonstrate  the  integrating  nature  of  the  EOSE  by  comparison  with  a  100  |im  Hertzian 
dipole  detector  that  has  been  shown  to  preserve  the  bipolar  nature  of  the  incident  electric  field 
waveform  [3].  Fig.  28a  shows  the  detected  waveform  with  the  dipole  detector  substituted  for  the 
EOSE  with  all  other  experimental  parameters  unchanged.  As  expected,  this  waveform  is  bipolar. 
The  time  integral  of  this  signal  is  shown  in  fig.  28b.  While  the  waveform  shapes  of  figs.  27  and 
28b  are  very  similar,  the  FWHM  time  of  the  former  signal  is  smaller.  To  understand  this,  we 
consider  the  response  time  of  the  EOSE,  neglecting  its  integrating  nature.  The  detection  response 
time  of  the  sampling  system  is  given  by  xs  =  ^Tpw  +  .  Here,  xpw  is  the  FWHM  optical 
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pulsewidth  of  the  probe  beam  and  is  equal  to  160  fs.  The  interaction  time  of  the  optical  pulse 
through  the  nonlinear  medium  Tit  =  2n£/c  is  110  fs  and  the  material  response  time  xm  is  <  10  fs 
[52].  Thus,  the  overall  FWHM  detection  response  time  ts  <  200  fs,  which  is  smaller  than  that  of 
our  dipole  detector. 

5.0  )  I  300 1 - 1 


Time  Delay  fps]  Time  Qe|ay  [pS] 

Figure  28a  Figure  28b 

It  should  be  possible  to  modify  the  sampling  element  so  that  the  incident  field  is  not 
diffracted.  This  can  be  accomplished  by  removing  the  electrodes  after  poling  or  by  using  an 
appropriate  organic  single  crystal.  In  addition  to  reducing  laser  noise,  improvement  in  the  detection 
sensitivity  should  be  possible  by  performing  the  experiment  in  a  transmission  geometry.  In  this 
geometry,  the  THz  field  would  again  be  normally  incident  on  the  EOSE.  The  probe  beam, 
however,  would  be  incident  at  the  phase  velocity  matching  angle.  This  angle  is  relatively  small  in 
polymers  due  to  the  low  refractive  index  dispersion  [44].  The  probe  beam,  therefore,  copropagates 
with  the  THz  beam  allowing  for  longer  interaction  lengths.  These  ideas  are  currently  being 
investigated. 

C.  Ultra-broadband  Coherent  Terahertz  Spectroscopy  Using  Optical  Rectification 
and  Electro-optic  Sampling: - 

In  recent  years,  remarkable  progress  has  been  made  in  the  development  of  spectroscopic 
capabilities  for  coherent  terahertz  (THz)  measurements.  A  key  ingredient  to  these  advances  is  the 
development  of  ultra-broadband  coherent  sources  and  detectors  of  THz  radiation.  Much  of  the 
work  reported  to  date  relies  on  the  use  of  photoconductive  elements  both  as  emitters  and  detectors. 
While  these  devices  are  highly  optimized  in  many  respects,  there  are  fundamental  limitations  in 
their  frequency  response  associated  with  the  natural  time  constants  for  carrier  dynamics.  The  ready 
availability  of  laser  pulses  with  durations  of  ~10  fs  suggests  the  potential  for  extending  the 
bandwidth  of  coherent  spectroscopy  to  significantly  higher  frequencies.  By  using  materials  with  an 
instantaneous  optical  response  for  both  emission  and  detection,  it  may  be  possible  to  capture  much 
of  this  enormous  bandwidth.  Thus,  the  use  of  materials  with  non-resonant  second-order  optical 
nonlinearities  is  optimally  suited  for  these  applications. 
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The  combination  of  broadband  optical  rectification  and  electro-optic  sampling  have  been 
demonstrated  in  %(2)  media  [53,  54].  In  order  to  enhance  the  sensitivity,  longer  interaction  lengths 
are  desirable,  necessitating  appropriate  consideration  of  phase  matching  constraints.  In  our  most 
recent  work,  we  have  addressed  this  issue  of  phase  matching  to  allow  for  greater  interaction 
between  the  optical  and  THz  pulses  in  the  nonlinear  medium  in  order  to  enhance  the  efficiency  of 
both  the  generation  and  detection  processes  [53],  We  have  demonstrated  that  dispersion  in  the 
optical  refractive  index  may  be  used  to  obtain  collinear,  non-critical  phase  matching  over  a  broad 
bandwidth  in  the  THz.  For  a  medium  with  dispersion  at  the  optical  frequencies,  the  phase  matching 
condition  for  both  difference-frequency  mixing  and  electro-optic  sampling  can  be  written  as 

„  r  (20) 

®THz 

This  relation  implies  that  phase  matching  is  achieved  when  the  phase  of  the  THz  wave  travels  at  the 
velocity  of  the  optical  pulse  envelope  (i.e.,  the  optical  group  velocity,  vg).  Optical  [55]  and  THz 
dispersion  [56]  data  indicate  that  sufficient  group  velocity  dispersion  exists  within  ZnTe  to  obtain 
phase  matching  at  fraz  ~  2  THz.  The  possibility  of  obtaining  phase  matching  in  this  manner  at  a 
pump  wavelength  of  800  nm  is  illustrated  in  fig.  29  for  ZnTe,  the  material  used  in  our  experimental 
investigations. 


Frequency  [THz] 


Figure  29 

The  coherence  decreases  because  of  dispersion  in  the  THz  refractive  index,  but  remains  large  over 
a  wide  range  of  frequencies.  Using  optical  rectification  and  electro-optic  sampling  in  <1 10>  ZnTe 
for  the  generation  and  coherent  detection  of  freely  propagating  THz  radiation,  respectively,  we 
have  demonstrated  spectral  sensitivity  beyond  3  THz  [53]. 

A  mode-locked  Ti:sapphire  laser  producing  130  fs  pulses  at  800  nm  with  a  76  MHz 
repetition  rate  was  used  to  generate  and  detect  the  FIR  pulses.  A  2  mm  diameter  pump  beam, 
appropriately  polarized  in  order  to  maximize  the  nonlinear  response,  was  chopped  and  used  to 
drive  a  0.9  mm  thick  <1 10>  ZnTe  emitter  at  normal  incidence.  The  THz  electromagnetic  transient 


T.F.  Heinz,  AFOSR  Final  Report 


45 


produced  by  optical  rectification  in  the  ZnTe  crystal  was  imaged  into  an  identical  ZnTe  crystal  (the 
E-0  sampling  element)  using  two  off-axis  paraboloidal  mirrors.  The  total  separation  between  the 
emitter  and  detector  was  50  cm.  A  pellicle  beamsplitter  was  interposed  in  the  THz  beam  line  to 
allow  for  co-propagation  of  the  optical  probe  and  THz  beams  through  the  EOSE.  The  temporal 
waveform  of  the  measured  THz  electric  field  is  shown  in  fig.  30.  The  full  width  at  half  maximum 
(FWHM)  pulsewidth  is  approximately  270  fs.  The  frequency  response  of  this  system  is  limited 
primarily  by  phase  matching  characteristics  (the  maximum  frequency  for  obtaining  a  coherence 
length  of  at  least  0.9  mm  is  ~2.5  THz).  In  addition  to  these  limitations,  the  dispersion  and 
absorption  caused  by  the  dominant  low  frequency  resonance  in  ZnTe  at  5.4  THz  strongly  attenuate 
the  high  frequency  components  [56]. 


Figure  30 

For  a  0.9  mm  thick  ZnTe  EOSE  and  shot-noise  limited  differential  detection, 
photomodulation  as  small  as  ~5  x  10"8  may  be  detected  (corresponding  to  a  minimum  detectable 
field  of  ~2  (mV/cm)/VHz).  Such  a  significant  improvement  in  sensitivity  of  the  EOSE  as 
compared  to  the  earlier  work  using  poled  polymers  [48]  is  largely  due  to  the  use  of  the  phase 
matched  transmission  geometry  which  allows  for  greater  interaction  lengths.  In  practice,  our 
sensitivity  is  -100  lower  and  is  believed  to  be  due  to  excess  laser  noise.  Modulation  of  the  pump 
beam  at  significantly  higher  frequencies  should  reduce  this  noise.  With  the  use  of  electro-optic 
materials  exhibiting  lower  dispersion  in  the  FIR,  it  should  be  possible  to  extend  the  bandwidth 
capability  of  this  system  with  enhanced  detection  sensitivity. 

VI.  Optoelectronic  Detection  of  Transient  Electric  Fields  in  High-Speed  Circuits 

The  rapid  progress  in  high  speed  electronic  devices  has  created  the  demand  for  new 
measurement  technologies.  For  maximum  utility,  these  techniques  must  exhibit  wide  bandwidth, 
accurate  signal  reproduction,  and  minimal  loading.  In  response  to  this  need,  several  ultrafast 
optoelectronic  methodologies  including  photoconductive  sampling,  electro-optic  sampling, 
photoemissive  sampling  ,  and  charge-sheet  probing  have  been  devised  [57].  Electro-optic  sampling 
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has  received  particular  attention  since  it  is  capable  of  providing  picosecond  temporal  resolution  and 
submillivolt  sensitivity  [46].  In  the  case  of  GaAs,  the  electro-optic  properties  of  the  substrate  itself 
may  be  exploited  to  perform  in-situ  measurements  [58].  However,  since  the  dominant  electronic 
material,  silicon,  does  not  exhibit  a  linear  electro-optic  effect,  an  external  probe  is  necessary. 

A.  High-Speed  Electrical  Sampling  Using  Optical  Second-Harmonic  Generation:- 

Optical  second  harmonic  generation  (SHG)  has  been  shown  to  be  an  effective  tool  for 
studying  surfaces  of  centrosymmetric  media  [59,  60].  A  high  degree  of  sensitivity  to  electric  fields 
has  been  demonstrated  in  studies  of  semiconductors,  both  with  fields  applied  normal  to  the  surface 
and  under  the  influence  of  naturally  occurring  depletion  fields  [61-63].  Recently,  Lupke  et  al. 
characterized  the  in-plane  electric-field  dependence  in  silicon  using  a  transmission  line  structure 

[64] ,  These  investigations  suggest  the  possibility  of  measuring  transient  electric  fields  in  high 
speed  silicon-based  devices.  We  have  recently  demonstrated  the  first  measurement  an  ultrafast 
electrical  pulse  propagating  on  a  transmission  line  with  subpicosecond  time  resolution  by  SHG 

[65] .  Since  the  method  exploits  the  nonlinear  optical  response  of  silicon  itself,  it  requires  no 
additional  probes  or  crystals  and  is  non-invasive.  In  this  approach,  we  expect  the  time  resolution  to 
be  limited  only  by  the  laser  pulse  duration  down  to  an  inherent  response  time  of  less  than  10  fs, 
while  the  spatial  resolution  is  limited  only  by  the  focusing  of  the  probe  beam.  Furthermore,  the 
present  scheme  is  readily  applicable  to  silicon-based  circuits  without  the  need  for  any  external 
probes  and  can,  therefore,  be  used  to  examine  high  density  integrated  circuits. 

We  fabricated  a  coplanar  waveguide  transmission  line  on  silicon-on-sapphire  that  was 
subsequently  ion  implanted.  In  addition  to  reducing  the  carrier  lifetime  to  ~0.6  ps,  we  have  found 
experimentally  that  this  0+  ion  implantation  step  significantly  improves  the  ohmic  nature  of  the 
contacts,  as  demonstrated  by  the  I-V  characteristics.  We  used  a  76  MHz  mode-locked  Ti:sapphire 
laser  operating  at  800  nm  with  150  fs  pulses  as  the  optical  source.  A  portion  of  the  laser  output 
served  to  launch  picosecond  electrical  pulses  onto  the  transmission  line  at  a  biased  photoconductive 
gap.  The  electrical  pulses  were  sampled  at  a  nearby  point  between  the  coplanar  lines  by  an  s- 
polarized  probe  beam  introduced  with  a  variable  temporal  delay.  For  the  SHG  measurement,  the  s- 
polarized  reflection  was  detected  after  passing  through  a  series  of  filters  using  a  photon  counter, 
since  this  configuration  suppressed  the  background  signal.  The  s-polarized  SH  intensity,  Is(2<0), 
may  be  related  to  the  DC  field,  Esdc,  oriented  along  the  s  direction  by  [60] 

lf“’“|3d2  +  5d2sEf|2(lf>)2,  (21) 

where  0)  is  the  frequency  of  the  fundamental  probe  beam  and  %(2)  and  %(3)  are  the  second  order  and 
third  order  susceptibilities,  respectively. 

First,  we  determined  the  dependence  of  the  SH  intensity  on  a  static  electric  field  produced 
by  applying  a  DC  bias  to  the  transmission  line  and  observed  a  parabolic  dependence  on  the  applied 
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voltage,  as  shown  in  fig.  31.  The  background  contribution  from  %(2)  is  seen  to  be  negligible,  as 
anticipated  for  ion  implanted  silicon. 


-10  -8  -6  -4  -2  0  2  4  6  8  10 

Bias  [volts] 

Figure  3 1 

We  then  measured  the  transient  electric  field  by  placing  a  bias  on  the  pump  finger  and  gating  the 
corresponding  photoconductive  gap.  The  time-resolved  SH  intensity  Is(2<0),  with  and  without  an 
applied  bias,  is  shown  in  fig.  32.  We  estimate  that  the  peak  pulse  amplitude  corresponds  to  a 
maximum  voltage  of  500  mV.  The  top  waveform  in  this  figure  has  a  time  response  (FWHM)  of 
approximately  1  ps.  As  expected,  the  SHG  derived  pulsewidth  is  comparable,  but  somewhat 
smaller,  than  the  electronic  cross-correlation  pulsewidth  of  1.55  ps  derived  from  two 
photoconductive  gaps. 


Time  Delay  [ps] 


Figure  32 

Electrical  pulses  propagating  on  such  coplanar  transmission  lines  have  been  characterized 
with  subpicosecond  time  resolution  and  100  mV/VHz  sensitivity.  Given  the  short  escape  depth 
(125  nm)  for  the  SH  signal  and  the  expected  electronic  nature  of  the  nonlinearity,  an  intrinsic 
response  time  of  <  10  fs  may  be  expected.  A  substantial  improvement  in  the  signal  to  noise  ratio 
may  be  expected  by  using  a  photon  energy  for  the  probe  lying  below  the  band  gap  of  Si.  In  this 
case,  the  technique  would  remain  non-perturbative  even  for  higher  probe  intensities,  since  no 
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carriers  will  be  produced  by  the  probe.  Since  the  observed  SH  signal  will  still  originate  from  the 
surface  region  of  the  silicon,  the  escape  depth  for  the  SH  radiation  (which  lies  above  the  band  gap) 
will  remain  short.  Furthermore,  significantly  tighter  focusing  of  the  incident  probe  beam,  coupled 
with  a  reduction  in  the  optical  pulse  duration,  will  yield  a  substantial  enhancement  in  the 
sensitivity. 

VII.  Interactions  and  Technology  Transfer 

During  the  course  of  the  research  described  in  this  report,  the  Columbia  group  has  had 
several  fruitful  interactions  and  collaborations  with  industrial  research  organizations,  government 
research  facilities  and  academic  institutions,  as  we  discuss  below. 

A.  Industrial  interactions 

The  Ultrafast  Optoelectronics  Group  at  Columbia  has  had  several  interactions  with  various 
researchers  at  AT&T  Bell  Laboratories.  The  ion-implanted  Si-on-sapphire  dipole  detectors  used  for 
our  coherent  THz  measurements  were  fabricated  with  the  assistance  of  Dr.  Peter  Smith  at  AT&T 
Bell  Labs,  Holmdel,  New  Jersey.  The  initial  work  on  using  large-aperture  Si  p-i-n  diodes  as  THz 
emitters  was  done  on  samples  provided  by  Dr.  Bahram  Jalali’s  group  at  AT&T  Bell  Labs,  Murray 
Hill,  New  Jersey. 

The  Columbia  group  has  collaborated  extensively  with  Dr.  James  T.  Yardley’s  organization 
at  the  AlliedSignal  Technology  Center  in  Morristown,  New  Jersey.  This  collaboration  has  allowed 
students  at  Columbia  to  have  full  access  to  materials,  processing  capabilities,  and  technical 
expertise  related  to  poled  polymers.  In  return,  this  interaction  has  yielded  valuable  information  to 
AlliedSignal  regarding  the  nonlinear  optical  properties  of  their  materials. 

The  Columbia  group  has  also  interacted  strongly  with  IBM  T.  J.  Watson  Research  Center 
in  Yorktown  Heights,  New  York.  In  response  to  an  interest  in  non-invasive  high-speed  probes  of 
circuits,  joint  work  was  initiated  with  Dr.  James  A.  Misewich  of  IBM  Yorktown  Heights.  These 
investigations  led  to  the  formulation  and  demonstration  of  the  SHG  technique  for  the  detection  of 
ultrafast  electrical  transients  that  is  described  in  the  previous  section  of  this  report.  Facilities  at  IBM 
Yorktown  Heights  have  been  made  available  for  fabrication  of  custom  optoelectronic  devices 
required  for  these  and  other  related  investigations.  Columbia  group  members  have  had  access  to 
IBM  facilities  and,  conversely,  IBM  staff  members  have  been  involved  in  ultrafast  optoelectronic 
measurements  carried  out  at  Columbia. 

B.  Interactions  with  DoD  labs  and  other  institutions 

The  Columbia  group  has  collaborated  extensively  with  Dr.  Jeffrey  D.  Morse  of  Lawrence 
Livermore  National  Labs  on  various  projects  involving  THz  emission  from  photoconducting 
devices.  Dr.  Morse’s  group  has  provided  various  custom-fabricated  semiconductor  samples  such 
as  low-temperature  MBE  grown  GaAs  and  GaAs  p-i-n  diodes. 
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With  regards  to  detailed  analyses  and  comprehensive  modeling  of  THz  emission  from  bulk 
semiconductors,  we  have  worked  closely  with  Dr.  A.V.  Kuznetsov  and  Prof.  C.J.  Stanton  from 
the  Department  of  Physics  at  the  University  of  Florida  at  Gainesville.  Our  group  has  also  interacted 
with  Dr.  L.  A.  Schlie  at  Phillips  Labs  in  Kirtland  Air  Force  Base,  and  has  performed  several 
measurements  using  available  femtosecond  laser  capabilities  to  characterize  the  operation  of  high¬ 
speed  THz  detectors  fabricated  at  Phillips  Labs. 
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